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Preface 


This  research  is  directed  at  improving  the  performance  of  regional  seismic 
discriminants  by  developing  new  approaches  for  reducing  the  scatter  in  discriminant 
measures  imparted  by  path-specific  waveguide  heterogeneity.  Conventional  procedures  for 
regionalization  of  discriminant  measures  involve  determining  empirical  distance- 
dependence  of  discriminants  for  regions  that  have  relatively  uniform  propagation 
characteristics,  with  no  blockage  of  phases.  Such  approaches  essentially  assume  that  the 
best  propagation  corrections  that  can  be  applied  are  those  for  a  single,  regionally  averaged 
empirical  distance  trend.  While  this  is  a  useful  and  necessary  step,  there  is  significant 
residual  variance  in  all  regional  discriminants  after  distance  correction,  and  this  scatter 
intrinsically  limits  the  performance  of  the  discriminants,  particularly  at  low  frequencies. 

Our  previous  work  has  demonstrated  that  consideration  of  individual  path 
properties  allows  further  empirical  reduction  of  discriminant  measure  variance,  using 
along-path  measures  obtained  from  data  bases  for  surface  topography,  crustal  thickness, 
and  sedimentary  layCT  thickness.  Several  path-specific  measures  are  found  to  have 
significant  correlations  with  regional  discriminant  measurements,  both  before  and  after 
correction  for  propagation  distance,  for  events  in  China  and  the  Middle  East.  For  station 
WMQ,  in  western  China,  variance  reductions  two  to  three  times  as  large  as  those  for 
distance  corrections  alone  can  be  achieved  by  simple  empirical  models  obtained  by 
multiple-regression  analysis.  We  are  applying  multiple-regression  analysis  for  several  path 
properties  to  all  CDSN  stations  and  all  broadband  stations  in  the  Middle  East  under  this 
contract  This  report  presents  two  preprints:  the  first  addresses  variance  reduction  of 
discriminant  measures  at  station  WMQ,  demonstrating  die  new  methodology,  and  the 
second  is  a  technical  study  comparing  discriminant  measurement  procedures  applied  to  data 
for  the  Middle  East  An  appendix  documents  the  full  set  of  empirical  relations  found  for 
our  data  for  WMQ. 
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ABSTRACT 


Short-period  regional  phases  play  an  important  role  in  identifying  low  magnitude 
seismic  events  in  the  context  of  monitoring  the  Comprehensive  Test  Ban  Treaty. 
Amplitude  ratios  of  regional  phases  comprised  mainly  of  P  wave  energy  (Pn,  Pg)  to  those 
comprised  mainly  of  S  wave  energy  (Sn,  Lg)  effectively  discriminate  between  explosions 
and  earthquakes  in  many  regions,  particularly  at  frequencies  higher  than  3  Hz.  At  lower 
frequencies,  discrimination  is  usually  poor  due  to  large  scatter  that  causes  overlapping  of 
event  populations.  Scatter  in  regional  discriminant  measures  such  as  Pg/Lg  ratios  is  caused 
by  both  source  and  propagation  effects,  and  reducing  the  scatter  imparted  by  the  latter  is 
essential  to  improving  the  discriminant  performance  when  events  do  not  share  identical 
paths.  Regional  phases  experience  distance-dependent  amplitude  variations  due  to  effects 
such  as  critical  angle  amplification,  geometric  spreading,  and  attenuation.  Discriminant 
measures  are  usually  corrected  for  empirically  determined  distance  trends  for  a  given 
region,  but  large  scatter  persists  after  such  corrections.  This  study  seeks  to  develop  more 
sophisticated  empirical  corrections  for  path  properties  in  order  to  further  reduce  the  scatter 
in  regional  discriminant  measures  caused  by  propagation  effects.  Broadband  seismic 
waveforms  recorded  at  station  WMQ,  in  western  China,  demonstrate  that  regional  phase 
ratios  such  as  Pn/Lg,  Pg/Lg,  Pn/Sn,  Pg/Sn,  and  Sn/Lg,  all  show  significant  distance 
dependence  for  frequencies  less  than  6  Hz.  However,  variations  in  crustal  structure  cause 
additional  path-specific  amplitude  fluctuations  that  are  not  accounted  for  by  regionally 
averaged  distance  corrections.  Blockage  of  Lg  phases  on  paths  traversing  the  margins  of 
the  Tibetan  Plateau  is  one  such  effect  Regression  analysis  demonstrates  that  regional 
discriminants  measured  at  WMQ  display  significant  correlations  with  path-specific 
properties  such  as  mean  elevation,  topographic  roughness,  basement  depth  and  crustal 
thickness.  Multiple  regressions  using  optimal  combinations  of  parameters  yield  corrections 
that  reduce  variance  in  Pg/Lg  measurements  for  frequencies  less  than  3  Hz  by  a  factor  of 
two  or  more  relative  to  standard  distance  corrections.  This  should  systematically  enhance 
the  performance  of  the  Pg/Lg  discriminant  at  low  firequencies.  The  method  presented  here 
can  be  used  for  all  regions  and  all  short-period  regional  discriminants.  It  is  likely  that  the 
extraordinary  crustal  heterogeneity  in  western  China  represents  an  extreme  case  of  path 
dependent  effects. 
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INTRODUCTION 


Seismic  monitoring  of  the  Comprehensive  Test  Ban  Treaty  (CTBT)  requires  high 
confidence  identification  of  sources  of  seismic  signals.  For  events  larger  than  about 
magnitude  4.0,  teleseismic  signals  usually  suffice  to  identify  the  source  (or  at  least  aUow 
ruling  it  out  as  a  possible  nuclear  explosion)  based  on  event  location,  source  depth,  or 
waveform  characteristics  such  as  P  wave  polarity  or  mj/Ms  ratios.  Events  with  lower 
magnitudes  usually  do  not  produce  adequate  teleseismic  signal  amplitudes  for  reliable 
source  identification,  particularly  at  longer  periods,  so  there  is  great  interest  in  exploiting 
the  larger  amplitude  short-period  signals  recorded  at  regional  distances  less  than  15(K)  km 
from  the  source.  Seismic  waves  at  regional  distances  tend  to  be  complicated  as  a  result  of 
multiple  reflections  and  phase  conversions  within  the  crustal  and  lithospheric  waveguide, 
particularly  in  the  short  period  range  (<  1  sec)  for  which  high  signal-to-noise  ratios  are 
observed.  The  primary  strategies  for  source  identification  with  regional  signals  involve 
characterirdng  frequency  dependence  of  either  particular  phases  (e.g.  high  frequency/low 
frequency  amplimdes  of  Pg  or  Lg)  or  of  the  overall  source  radiation  (e.g.  Ml/M„  ratios),  or 
characteri2dng  the  relative  amounts  of  S  and  P  energy  released  by  the  source  (e.g.  Pg/Lg 
ratios).  The  latter  approach  is  particularly  promising  as  it  can  exploit  relatively  narrow 
band,  high  frequency  observations  which  have  good  signal-to-noise  properties,  and  the  P/S 
ratios  reflect  source  radiation  effects  that  directly  distinguish  earthquakes  from  explosions. 

There  are  four  major  regional  phases  that  are  commonly  considered  for  regional  P/S 
measurements;  Pn  and  Sn,  which  are  head  waves  and  lithospheric  guided  waves  that 
propagate  primarily  in  the  upper  mantle  lid  just  below  the  Moho  boundary,  and  Pg  and  Lg, 
which  involve  direct  P  or  S  crustal  arrivals  as  well  as  an  extended  coda  of  internal  crustal 
reflections  and  reverberations.  PnA^g,  PgA^g,  Pn/Sn,  and  Pg/Sn  ratios  have  aU  been 
considered  for  explosion  discrimination,  with  one  ratio  or  another  usually  being  favored 
from  region  to  region  depending  on  local  crustal  properties  that  influence  energy 
partitioning  in  the  waveguide  (some  regions  excite  very  weak  Pn  or  Sn  arrivals,  while  m 
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other  areas  these  may  be  stronger  than  Pg  or  Lg).  Experience  has  shown  that  these  P/S 
ratios  tend  to  give  poor  discrimination  between  explosion  and  earthquake  signals  for 
frequencies  near  1  Hz  when  standard  distance  coirections  are  applied,  but  as  the  frequency 
increases  the  separation  of  explosion  and  earthquake  populations  improves  rapidly,  often 
proving  quite  diagnostic  for  frequencies  above  3  to  5  Hz  (e.g.  Waller  et  al.,  1995;  Taylor, 
1996;  Hartse  et  al.,  1996).  Regional  P/S  type  ratios  have  significant  scatter  for  both 
earthquake  and  explosion  populations  even  on  a  common  path  (e.g.  Walter  et  al.,  1995), 
which  indicates  that  source  effects  such  as  depth,  focal  mechanism,  stress  drop,  and 
proximity  to  the  water  table  do  influence  the  measures.  This  scatter  is  difficult  to  reduce 
other  than  by  averaging  multiple  observations  for  a  given  event,  which  can  reduce  effects 
such  as  those  caused  by  radiation  patterns.  However,  the  scatter  in  regional  phase  ratios  is 
significantly  greater  when  diverse  propagation  paths  are  considered  together.  This  is  the 
typical  case  for  CTBT  monitoring  situations,  in  which  sparse  numbers  of  stations  record 
events  spread  over  large  areas  at  regional  distances.  Accounting  for  propagation  effects  on 
regional  discriminant  measurements  is  necessary  if  the  scatter  is  to  be  reduced  to  the  level 
of  intrinsic  source  variability,  which  is  the  best  we  can  hope  for.  Low  false-alarm 
identification  of  outlier  signals  that  may  represent  violations  of  the  CTBT  requires  tightly 
defined  discriminant  populations  for  observations  on  many  different  paths. 

Efforts  to  understand  wave  propagation  effects  on  regional  phases  involve  a  wide 
variety  of  observational  and  theoretical  studies,  but  we  are  still  far  from  having  good 
predictive  capabilities  for  most  regional  signals.  This  is  because  the  crustal  and  lithospheric 
waveguide  is  very  complex  with  multiple  scales  of  heterogeneity.  Our  most  detailed  crustal 
models  are  still  very  simple  and  poorly  constrained,  and  many  areas  of  CTBT  monitoring 
concern  having  totally  unknown  crustal  structures,  and  high  frequency  wave  propagation  in 
the  three-dimensional  heterogeneous  crust  is  poorly  understood  and  cannot  be  accurately 
modeled  at  present  for  useful  propagation  distances.  As  a  result,  deterministic  approaches 
to  calibrating  regional  discriminants  and  conecting  observations  for  propagation  effects  are 
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very  limited,  mainly  being  restricted  to  source  modeling  of  longer  period  energy  (>3  sec)  in 
very  well-studied  regions  and  for  quantifying  effects  such  as  phase  blockage  by  waveguide 
disruptions. 

Most  efforts  to  develop  regional  discriminants  rely  on  empirical  approaches  for 
reducing  propagation  effects,  following  the  basic  strategy  that  has  long  been  used  for 
seismic  magnitudes.  This  involves  determining  empirical  curves  of  amplitude  observations 
as  a  function  of  epicentral  distance,  which  are  used  to  equalne  individual  phase  amplitudes 
or  amplitude  ratios  to  a  common  distance,  “correcting”  for  propagation  effects.  Typically, 
the  distance  corrections  are  determined  for  a  well-defined  region,  recogmzing  that  crustal 
variations  strongly  affect  the  amplitudes  of  regional  phases.  Amplitude-distance  or 
amplitude  ratio-distance  relationships  defined  for  a  given  region  in  this  fashion  implicitly 
assume  that  propagation  effects  (or  differential  propagation  effects  in  the  case  of  ratios)  can 
be  characterized  by  a  single  function  of  path  length  alone.  This  is  the  same  as  asserting  that 
the  crust  throughout  the  region  is  laterally  uniform,  such  that  all  paths  have  the  same 
behavior  as  long  as  their  path  length  is  the  same.  Unfortunately,  for  regional  phases 
(unlike  teleseismic  magnitudes)  this  is  a  poor  assumption  in  almost  all  regions  of  the  world, 
especially  in  tectonically  active  areas  where  natural  seismicity  presents  great  challenges  for 
CTBT  monitoring. 

There  is  certainly  value  in  determining  a  regionally  averaged  distance  dependence 
for  a  given  discriminant  measure.  Differential  geometric  spreading,  attenuation  and  average 
crustal  waveguide  effects  on  seismic  energy  partitioning  can  be  accounted  for  to  first  order. 
This  assertion  is  borne  out  by  practice:  almost  all  discriminant  measures  exhibit  some 
empirical  distance  dependence,  at  least  in  lower  frequency  bands.  However,  it  is  clearly 
desirable  to  account  for  more  realistic  propagation  effects  by  considering  specific  path 
properties,  which  cause  fluctuations  about  the  average  pattern  mamfested  in  regional 
distance  corrections.  Lacking  deterministic  capabilities  for  this,  and  because  the  problem 
does  not  reduce  to  a  simple  station  correction  t3qie  of  procedure  as  for  teleseismic 
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magnitudes,  we  pursue  empirical  approaches  that  have  been  motivated  by  many  previous 
studies  that  document  the  importance  of  specific  path  properties. 

For  example,  it  has  been  demonstrated  that  laterally  heterogeneous  velocity 
structure,  and  changes  in  crustal  discontinuities  and  crustal  thickness  affect  regional  phases 
(e.g.  Kennett,  1986,  1989;  Zhang  et  al,  1994).  Blockage  of  Lg  phases  by  segments  of 
oceanic  crust  as  short  as  100  to  200  km  is  one  classic  example  (e.g.  Press  and  Ewing, 
1952;  Zhang  and  Lay,  1995).  Lg  blockage  has  also  been  noted  in  continental  areas  with 
large  changes  in  crustal  thickness,  as  for  the  Tibetan  Plateau  (Ruzaikin  et  al.,  1911  \  Ni  and 
Barazangi,  1983;  McNamara  et  al.,  1996)  and  near  the  Alpine  Mountains  (Campillo  et  al., 
1993).  Sn  phases  are  strongly  attenuated  in  localized  areas  of  the  northern  Tibetan  Plateau 
(Ni  and  Barazangi,  1983)  and  in  the  northwestemmost  Iranian  and  Turkish  Plateaus  and 
between  the  Black  Sea  and  Caspian  Sea  (Kadinsky-Cade  et  al.,  1981;  Rodgers  et  al., 
1996).  Structure  on  the  Moho  has  been  foimd  to  affect  Pn  phases  in  Scandinavia  (Kvaema 
and  Doombos,  1991).  Topographic  characteristics  are  sometimes  correlated  with  regional 
phase  behavior  as  well  (Baumgardt,  1990;  Zhang  and  Lay,  1994;  Zhang  et  al.,  1996),  as 
are  sedimentary  structures  (Baumgardt  and  Der,  1994;  Zhang  et  al.,  1994).  While  it  is 
possible  in  practice  to  identify  regions  of  phase  blockage  and  to  apply  alternate 
discrimination  procedures  fliat  are  not  biased  by  such  effects  in  those  areas,  the  many 
studies  mentioned  above  establish  that  waveguide  irregularities  are  expected  to  be 
ubiquitous  and  resulting  fluctuations  in  discriminant  measures  will  exist  in  every  region. 
Even  effects  such  as  blockage  may  be  cumulative  effects,  with  a  continuum  of  behavior  that 
can  not  be  readily  regionalized.  Empirically  reducing  the  effects  seems  to  be  the  best 
approach,  much  as  is  done  for  standard  distance  corrections. 

In  this  study,  we  pursue  the  empirical  calibration  approach  suggested  by  Zhang  and 
Lay  (1994)  and  Zhang  etal.  (1994, 1996).  These  studies  used  the  available  data  bases  on 
crustal  properties  in  different  regions  to  develop  path-specific  parameters  which  are  tested 
for  correlation  with  regional  phase  amplitude  ratios.  The  most  reliable  parameters  that  can 
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be  deteraiined  for  any  path  are  those  involving  surface  topography,  for  which  highly 
precise  models  are  available.  Crude  models  of  crustal  thickness  and  sedimentary  layer 
thickness  are  also  available  for  most  of  Eurasia.  These  models  are  used  parametrically  as 
surrogates  for  the  complex  propagation  effects  produced  by  the  heterogeneous  waveguide. 
Ratios  of  regional  phases  recorded  at  station  WMQ,  in  western  China,  are  analyzed  for 
regionally-averaged  distance  effects  as  well  as  path-specific  effects  parameterized  by 
measures  of  surface  topography,  sedimentary  layer  thickness  and  crustal  thickness.  Both 
point-wise  properties  (e.g.,  maximum  or  minimum  crustal  thickness  on  each  path)  and 
path-integrated  prc^erties  (e.g.,  mean  altitude,  mean  crustal  thickness)  are  considered.  Our 
goal  is  to  reduce  scatter  in  fiequency  dependent  P/S  type  seismic  discriminants  by 
performing  multiple  regressions  to  establish  empirical  models  for  propagation  corrections. 
Results  are  presented  for  Pg/Lg  observations.  The  models  are  not  unique,  and  they  do  not 
reveal  the  causal  physics  controlling  the  regional  phases,  but  the  same  statements  are  true  of 
empirical  distance  corrections  (which  are  based  on,  at  most,  very  simplified  models  of 
geometric  spreading  and  attenuation  that  are  usually  not  even  consistent  with  our  gross 
knowledge  of  Earth  layering  and  regional  phase  propagation).  Significant  variance 
reductions  can  be  achieved  with  plausible  combinations  of  waveguide  parameters,  enabling 
improved  discriminant  performance  at  lower  frequencies.  Such  procedures  can  be  readily 
incorporated  into  CIBT  monitoring  operations.  The  empirical  relations  that  are  obtained 
also  provide  a  starting  point  for  quantitative  modeling  efforts  that  can  be  pursued  to 
quantify  the  underlying  wave  propagation  effects. 

DATA 

Western  China  is  a  region  with  very  diffuse  seismicity  (Figure  1),  sparse  seismic 
instrumentation,  acute  crustal  heterogeneity,  and  prior  nuclear  testing  experience.  Thus,  it 
is  a  region  of  interest  for  CTBT  monitoring,  and  many  of  the  practical  issues  about 
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Figure  1.  Seismicity  in  western  China  and  its  vicinity.  The  map  shows  the  epicentral 
locations  for  earthquakes  from  1988  to  1995  from  the  PDE  catalogue.  The  small  dots 
represent  epicenters  of  earthquakes  with  magnitudes  4.0  <  mb  <  5.0;  the  small  open  circles 
are  epicenters  of  earthquakes  with  magnitudes  of  5.0  ^  mb  <  6.0;  and  large  open  circles  are 
epicenters  for  earthquakes  with  magnitudes  of  6.0  <  mb. 


calibration  of  regional  discriminants  are  relevant  to  this  region.  We  analyze  87  earthquakes 
with  magnitudes  of  4.5  <  mb  ^  6.1  that  occurred  between  1988  and  1995  in  western  China 

and  its  vicinity.  These  events  were  all  recorded  on  the  broadband  vertical  component  at 
station  WMQ,  with  the  path  coverage  being  shown  in  Figure  2.  Figure  3  shows  a  few 
example  seismograms,  indicating  the  great  variability  in  broadband  waveform  complexity 
on  the  diffaent  paths.  To  avoid  ambiguity  in  regional  phase  recognition  at  very  close 
distances,  only  events  with  epicentral  distances  greater  than  240  km  are  used,  with  the  most 
distant  observations  being  2100  km  from  the  station.  The  events  are  aU  large  enough  to 
have  Harvard  centroid  moment  tensor  (CMT)  focal  mechanisms,  and  they  have  fairly  well 
determined  locations  as  a  result  of  numerous  teleseismic  detections.  For  events  prior  to 
1994  we  use  locations  given  by  the  International  Seismological  Centre  (ISC)  bulletins,  and 
for  events  after  1994  we  use  locations  from  the  USGS  Preliminary  Determination  of 
Epicenters  (PDE)  catalog.  Only  events  with  focal  depths  less  than  50  km  are  used  in  this 
study,  but  there  are  relatively  large  uncertainties  in  source  depth  due  to  the  sparse  regional 
station  coverage  and  the  complex  crustal  and  upper  mantle  structure  in  the  region. 

Following  Walter  era/.  (1995)  and  Hartse  etal.  (1996),  measurements  of  regional 
seismic  phases  are  made  for  all  ground  motions  in  specific  windows.  The  Pn  window 

begins  at  the  onset  of  the  signal  (aH  events  are  at  distances  beyond  the  crossover  distance) 
and  ends  at  the  start  of  the  Pg  window.  The  Pg  and  Sn  windows  are  between  group 
velocities  of  6.2  to  5.2  km/sec  and  4.8  to  4.0  km/sec,  respectively.  The  Lg  window 
extends  from  a  group  velocity  of  3.6  km/sec  to  a  group  velocity  of  3.0  km/sec.  We 
carefully  checked  the  validity  of  these  windows  for  capturing  the  primary  energy  in  each 
signal.  For  each  phase,  the  vertical  component  broadband  records  (approximately 
proportional  to  ground  velocity)  were  windowed  and  band  pass  filtered  into  four  frequency 
bands:  0.75  to  1.5  Hz,  1.5  to  3.0  Hz,  3.0  to  6.0  Hz,  and  6.0  to  9.0  Hz,  using  Butterworth 
filters.  We  did  not  deconvolve  the  instrument  response,  as  we  consider  only  ratios  of 
phases  in  the  same  frequency  band,  which  explicitly  cancels  the  instrument  effect. 
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Seismic  Waveforms  Recorded  at  WMQ 


0  200  400  600 


TIME  (s) 

Figure  3.  Examples  of  broadband  vertical  component  seismograms  recorded  at  station 
WMQ  for  events  in  different  source  regions.  The  epicentral  distance  from  WMQ  is  given 
by  Dist,  and  BAZ  indicates  the  azimuth  from  the  station  toward  the  event. 


Previous  studies  have  demonstrated  that  rms  measures  of  regional  phases  like  Lg  have 
remarkable  stability  due  to  the  averaging  of  multiple  arrivals  (Ringdal  et  al.  1992;  Gupta  et 
al.,  1992),  so  we  compute  rms  amplitudes  of  each  phase  in  the  corresponding  filtered 
group  velocity  windows.  Rodgers  et  al.  (1997)  have  considered  the  various  regional  phase 
measurement  procedures,  and  while  they  recommend  frequency  domain  spectral  averaging, 
the  results  of  time  domain  measures  of  velocity  traces  are  not  significantly  different  than  for 
most  other  commonly  used  measures  (rms  ground  displacements,  spectral  averages,  etc.). 
Noise  levels  are  computed  using  rms  measures  of  the  filtered  seismograms  in  the  10  sec 
window  preceding  the  manually  picked  first  arrival  time.  Events  with  Pn  signal-to-noise 
ratios  greater  than  two  are  considered  to  be  acceptable,  given  that  Pn  is  usually  the  weakest 
regional  phase  in  these  data,  except  when  Lg  blockage  occurs.  For  most  of  the  events,  die 
signal-to-noise  ratios  are  between  10  to  100. 

The  ray  paths  coimecting  the  earthquakes  and  WMQ  sample  a  large  area  in  the 
vicinity  of  western  China,  including  some  of  the  most  dramatic  topographic  features  in  the 
world  such  as  the  Tibetan  Plateau,  and  the  Tien  Shan,  Pamir  and  Hindu  Kush  Mountains. 
Many  ray  paths  project  along  or  across  the  mountain  belts  or  the  margins  of  the  Tibetan 
Plateau.  Topography  in  this  region  is  illustrated  in  Figure  4.  The  topographic  relief  data 
are  from  the  GLOBE  (Global  Land  One-km  Base  Elevation)  project  with  a  spatial 
resolution  of  30  arc-second  (one-km).  This  database  was  developed  by  the  U.S.  Defense 
Mapping  Agency  (DMA)  and  made  available  by  the  National  Geophysical  Data  Center 
(NGDC)  of  NOAA.  The  lowest  elevation  is  -152  m  near  the  Tinfan  Basin  (42.6°N, 
89.3°E).  The  highest  elevation  occurs  in  the  Himalayas,  where  many  mountain  peaks  are 
well  above  8000  m  in  altitude.  Surface  topography  is  the  most  precisely  known  attribute  of 
the  regional  crustal  structure,  thus  using  it  as  a  guide  to  wave  propagation  influences  is 
clearly  desirable.  The  general  notion  is  that  topography  reflects  the  tectonic  character  of  Ihe 
crust,  with  rough  topography  in  active  regions  with  disrupted  crust,  high  topography  in 
areas  of  tectonically  thickened  crust,  laterally  variable  topography  in  areas  of  rapidly 
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changing  crust  and  mantle  structure,  and  so  on.  Given  variable  degree  of  isostatic 
compensation  and  smoothing  effects  of  sedimentation,  topography  cannot  be  expected  to 
reflect  all  crustal  structure  variations,  but  it  is  clearly  a  starting  point. 

Much  less  is  known  about  the  internal  crustal  structure  of  the  Eurasian  crust,  but 
low  resolution  models  of  crustal  thickness  (Moho  depth),  as  shown  in  Figure  5,  and 
sediment  thickness,  as  shown  in  Figure  6,  are  available  from  the  database  compiled  by 
Cornell  University.  These  data  were  digitized  from  maps  compiled  by  Kunin  and  Sheykh- 
Zade  (1983)  and  Kunin  (1987)  based  on  deep  seismic  sounding  (DSS)  profiles  and  gridded 
on  a  10  km  grid  (Fielding  et  al.,  1992).  There  is  no  question  that  the  accuracy  is  limited, 
but  the  maps  were  derived  independently  from  surface  topography,  and  when  averaged 
along  any  given  raypath  they  probably  provide  first  order  average  path  properties  that  are 
not  greatly  in  error.  There  is  much  less  confidence  in  any  point-wise  property  or  local 
gradient  in  the  properties  in  these  models,  and  we  keep  this  in  mind  when  using  these 
characterizations  of  the  waveguide. 

Profiles  of  topography,  Moho  depth,  and  sedimentary  layer  thickness  along  each 
great  circle  path  from  event  epicenter  to  station  WMQ  were  extracted  using  software  from 
Wessel  and  Smith  (1991).  In  order  to  reduce  these  functions  to  measures  that  can  be 
correlated  with  regional  phase  amplitude  ratios,  we  must  compute  summary  parameters. 
For  topography,  which  is  very  precisely  known,  we  computer  four  path-specific 
parameters:  mean  elevation,  rms  roughness  (variance)  of  topography,  rms  slope  of 
topography  from  point  to  point  along  the  path,  and  topographic  skewness,  or  the  third 
moment  of  topography,  which  characterizes  the  degree  of  asymmetry  of  topography  around 
its  mean  value  (Bennett  and  Mattsson,  1989;  Zhang  and  Lay,  1994).  We  also  calculate 
minimum,  maximum,  and  mean  crustal  thickness,  and  maximum  and  mean  sediment 
thickness  along  each  path.  The  minimum  and  maximum  measures  are  local  values,  and 
hence  subject  to  much  greater  uncertainties  than  the  path-averaged  mean  values,  but  local 
extremes  may  be  responsible  for  important  effects  such  as  blockage  of  phases.  Other 
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Figure  5.  Depth  to  the  Moho  discontinuity  in  the  study  area.  HK,  Hindu  Kush 
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Figure  6  S^mentaiy  layer  thickness  in  the  study  area.  QB,  the  Qsaidam  Basin;  JB,  the 

Junggar  Basin;  SB,  the  Sichuan  Basin. 


databases  could  be  considered,  such  as  path  averages  through  tomographic  models  derived 
from  surface  waves  or  attenuation  models,  but  at  present  these  are  all  even  less  reliable  for 
western  China  than  the  gross  crustal  models,  so  we  restrict  our  attention  to  the  measures 

listed  above. 

PATH  EFFECTS  FOR  REGIONAL  DISCRIMINANTS 

We  have  systematically  analyzed  the  path  effects  on  all  Pn/Cg,  Pg/Lg>  Pn/Sn> 
Pg/Sn,  and  Sn/Lg  ratios  measured  in  different  frequency  bands  for  the  87  events  recorded 
at  WMQ.  The  first  effect  to  be  considered  is  simple  propagation  path  length.  The  relative 
amplitudes  of  regional  phases  experience  a  general  distance-dependent  variation  caused  by 
differential  attenuation  and  geometric  spreading  (Kennett,  1989;  McCormack  et  al.,  1995; 
Zhang  et  al.,  1996).  This  distance  dependence  is  expected  to  be  very  complex,  if  one  uses 
regional  phase  calculations  for  a  simple  crustal  model,  but  such  calculations  are  usually 
poor  approximations  to  actual  data,  which  have  many  extra  airivals  caused  by  scattering 
and  complex  crustal  layeiing.  In  practice  it  is  found  that  the  regional  phase  signals, 
especially  when  measured  as  rms  values  over  a  group  velocity  window  which  contains 
multiple  arrivals,  usually  have  overall  amplitude  patterns  that  can  be  pretty  well  represented 
by  a  power  law  (A  form  (Sereno,  1991)  or  an  exponential  form  (e<  ^'^^)  (Baumgardt  and 
Der,  1994).  Zhang  and  Lay  (1996)  proposed  a  slightly  different  exponential  form  for 
characterizing  distance  dependence  of  phase  ratios,  R(A)=10^^,  which  is  adopted  here. 
There  is  not  a  great  deal  of  difference  in  these  parameterizations,  and  similar  results  are 
obtained  even  if  the  corrections  are  based  on  log(A)  rather  than  A,  as  in  some  studies.  It  is 
important  to  note  that  even  if  the  crust  were  laterally  uniform,  these  simple  distance 
expressions  are  not  strictly  correct,  as  the  nature  of  most  regional  phases  is  that  their 
amplitudes  experience  complex  variations  with  distance  due  to  critical  angle  effects  (for 
example,  multiple  SmS  arrivals  in  Lg  phases  can  produce  a  ‘bursty’  pattern  with  distance). 
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Nonetheless,  simple  amplitude-distance  parameterizations  like  these  are  the  norm  and  they 
work  very  well,  as  shown  below.  Given  our  limited  understanding  of  crustal  structure  in 
most  regions,  using  these  simple  forms  is  defensible,  albeit  limited  to  a  veiy  gross 
characterization  of  propagation  effects. 

Linear  regressions  for  all  of  the  regional  phase  ratios  measured  at  WMQ  show  clear 
distance  dependence.  Apparently,  when  the  rms  amplitude  averaging  of  each  phase 
window  is  performed,  the  progressively  accumulating  relative  attenuation  and  geometric 
spreading  effects  prove  to  be  primary  factors  controlling  regional  phase  energy  ratios.  The 
highest  correlations  with  propagation  distance  (up  to  0.87)  are  found  for  Log(Pn/Sn),  with 

the  frequency  dependent  measures  being  shown  in  Figure  7.  This  is  likely  due  to  the  fact 
that  both  Pn  and  Sn  phases  mainly  propagate  within  the  upper  mantle  lid  beneath  the  Moho, 

and  relative  attenuation  and  geometric  spreading  differences  associated  with  varying  P  and 
S  velocity  gradients  in  the  lid  produce  smoothly  accumulating  effects  with  distance 
traveled.  The  slopes  of  the  linear  regressions  and  the  linear  correlation  coefficients 
diminish  with  increasing  frequency,  as  proves  true  for  all  amplitude  ratios.  Log(Pn/Lg) 

ratios  also  have  strong  correlations  with  distance  for  frequencies  up  to  3  Hz  (>0.80;  Figure 
A3a).  Other  amplitude  ratios  of  Pg/Lg,  Pg/Sn,  and  Sn/Lg  also  display  significant  distance 

dependence  with  the  correlation  coefficients  ranging  between  0.4  to  0.7  in  the  three 
passbands  below  6  Hz  (see  Appendix  Figures  A2a,  A4a  and  A6a). 

While  the  amplitude  ratios  of  various  regional  phases  show  strong  correlations  with 
path  length,  independent  relationships  are  also  found  with  topographic  parameters  and 
other  waveguide  characteristics.  For  the  various  topographic  parameters  considered,  mean 
path  elevation  and  average  surface  roughness  along  each  path  have  the  most  systematic 
correlations,  as  has  been  observed  in  other  regions  (Zhang  and  Lay,  1994;  Zhang  et  al., 
1994, 1996).  The  correlations  for  the  5  frequency  dependent  amplitude  ratio  measurements 
with  these  parameters,  without  correction  for  path  length,  are  shown  in  Figures  A2b,  A2d, 
A3b,  A3d,  A4b,  A4d,  A5b,  A5d,  A6b  and  A6d.  In  every  case  the  linear  correlation 
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Figure  7.  Log(Pn/Sn)  as  a  function  of  distance  from  the  station  in  various  frequency 
bands.  CC,  linear  correlation  coefficients;  SD,  standard  deviation  of  the  linear  regression; 
K,  slope  of  the  regression  line.  The  dotted  lines  indicate  ±  1  standard  deviation  about  the 
regression  line. 
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coefficients  are  larger  than  0.5  for  frequency  bands  0.75-1.5  Hz  and  1.5-3  Hz,  and  below 
0.5  for  higher  frequencies.  In  fact,  for  all  ratios  other  than  those  involving  Pn,  the  linear 
correlation  coefficients  with  mean  path  elevation  and  surface  roughness  are  larger  than 
those  with  path  length  for  the  two  lower  frequency  bands.  The  path-averaged  rms  slope  of 
topography  shows  little  correlation  (linear  correlation  coefficients  less  than  0.2)  with  the 
logarithmic  amplitude  ratios  in  all  cases  (Figures  A2j,  A3j,  A4j,  A5j,  A6j),  unlike  the 
results  for  the  Russian  platform  reported  in  Zhang  and  Lay  (1994).  A  much  larger  range  of 
rms  slope  variations  is  spanned  by  our  data  (from  20  to  220  m  per  km)  in  contrast  to  the 
small  range  (2  to  20  m  per  km)  in  the  study  of  Zhang  and  Lay  (1994),  so  their  result 
appears  to  be  an  isolated  case.  The  fourth  topography  parameter,  skewness,  shows 
moderate  negative  correlations  (linear  correlation  coefficients  of  -0.33  to  -0.56)  with  tire 
regional  phase  amplitude  ratios  for  frequencies  below  3  Hz  (Figures  A2h,  A3h,  A4h,  A5h 
and  A6h). 

Zhang et  al.  (1996)  found  that  in  the  western  U.S.,  the  strongest  single-parameter 
empirical  correlations  were  found  for  the  product  of  mean  elevation  times  path  length  or 
mean  roughness  time  path  length,  indicative  of  a  cumulative  effect  on  the  amplitude  ratios. 
This  is  also  true  of  the  WMQ  observations,  with  the  linear  correlation  coefficients  being 
systematically  higher  for  frequencies  lower  than  3  Hz  for  Pg/Lg  (Figures  A2f,  A2g)  Pn/Lg 
(Figures  A3f,  A3g),  Pg/Sn  (Figures  A4f,  A4g),  and  Sn/Lg  (Figure  A6f,  A6g).  One 
exception  is  for  Log(Pn/Sn),  for  which  the  parameter  distance  times  roughness  does  not 

improve  the  correlations  relative  to  distance  alone.  This  suggests  that  topographic 
roughness  has  litfle  affect  on  the  Moho  head  waves,  as  might  be  expected.  For  the  two 
lower  fiequency  bands  of  Pg/Lg  measurements  the  linear  correlation  coefficients  are  higher 
than  0.79  for  path  parameters  based  on  the  products  of  distance  and  mean  elevation  or 
distance  and  surface  roughness.  This  can  be  compared  to  distance  correlations  of  0.56  to 
0.66  to  infer  that  there  are  significant  variance  reductions  to  be  achieved  by  including 
additional  path  information.  The  most  dramatic  correlations  are  for  Pn/Lg  measurements 
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with  these  same  products  of  distance  times  mean  elevation  or  roughness,  with  linear 
correlation  coefficients  greater  than  0.88,  although  for  Pn/Lg  the  distance  trends  alone  are 

very  strong  (linear  correlations  larger  than  0.8). 

Proceeding  to  additional  path  characteristics,  we  anticipate  that  crustal  thickness 
may  be  an  important  parameter  because  it  provides  a  direct  measure  of  the  overall  shape  of 
the  waveguide  in  which  the  regional  seismic  energy  is  partitioned.  The  WMQ  data  are 
valuable  for  assessing  this  given  the  large  variations  in  crustal  thickness  in  the  region 
(Figure  5).  Seismic  reflection  and  refraction  studies  have  demonstrated  that  crustal 

thickness  in  the  Tibetan  Plateau  is  about  65  to  75  km  (Him  et  al.,  1984),  making  it  one  of 
the  thickest  cmstal  areas  in  the  world.  Lg  waves  are  apparently  highly  attenuated  due  to  a 
low  Lg  Q  value  within  the  Tibetan  Plateau  and  blocked  at  its  northern  and  southern 

boundaries  because  of  changes  in  cmstal  thickness  across  the  margin  (McNamara  et  al, 
1996).  We  find  that  correlations  of  Log(Pg/Lg)  with  mean  (Figure  8)  and  maximum 

(Figure  A2o)  cmstal  thickness  are  higher  than  those  with  distance  or  topography  for  tiie 
0.75-1.5  Hz  band,  and  comparable  for  the  1.5-3  Hz  band.  The  correlations  are  very  low 
for  higher  frequencies.  To  a  large  extent,  the  low  frequency  trends  reflect  the  partial  and 
complete  blockage  of  Lg  phases  traversing  the  margin  of  the  Tibetan  Plateau,  which  is  also 
partly  responsible  for  the  trends  with  topography  (note  the  strong  correlations  between 
Figures  4  and  5).  This  is  consistent  with  modeling  efforts  that  used  ray  diagrams  and 
topography  in  central  Asia,  including  Tibet,  to  demonstrate  that  Pg  and  Lg  propagation  are 
affected  by  changes  in  the  cmstal  thickness  (Kennett,  1989;  Bostock  and  Kennett,  1990). 
The  correlations  appear  to  be  stronger  at  lower  fi:equencies  because  low  frequency  Lg  has 
more  predictable  effects  than  the  strongly  scattered  high  frequency  energy  (e.g.,  Shih  et  al, 
1994).  The  correlation  coefficients  for  Log(Pg/Sn)  uud  mean  cmstal  thickness  (Figure 

A41)  are  slightly  higher  than  those  with  distance  (Figure  Ana),  but  slightly  lower  than  those 
with  mean  elevation.  Correlations  of  mean  or  maximum  cmstal  thickness  with  Pn/Lg, 
Pn/Sn,  and  Sn/Lg  are  comparable  to,  or  lower  than  those  with  topography  and,  in  general. 
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Figure  8.  Log(Pg/Lg)  as  a  function  of  mean  crustal  thickness  in  various  frequency  bands. 
CC,  linear  correlation  coefficient;  SD,  standard  deviation  of  the  linear  regression;  K,  slope 
of  the  regression  line.  The  dotted  lines  indicate  ±  1  standard  deviation  about  the  regression 


weaker  than  those  with  distance  (Figures  A31,  A3o,  A51,  A5o,  A61,  A6o).  In  every  case 
we  found  weak  correlations  (<0.45)  with  minimum  crustal  thickness  on  the  path  (Figures 
A2n,  A3n,  A4n,  A5n,  A6n)  so  this  is  not  a  very  effective  parameter  for  characterizing 

regional  phase  amplitude  variations  in  this  region. 

Sedimentary  structures  may  cause  blockage  of  Lg  by  capturing  energy  or  by 
scattering  it  out  of  the  waveguide  (e.g.  Baumgardt  and  Der,  1994;  Zhang  et  aL,  1994),  and 
some  of  our  paths  do  traverse  areas  of  thick  sediments  in  the  Tarim  Basin,  the  Junggar 
Basin,  the  Sichuan  Basin  ,  the  Qsaidam  Basin,  and  Qrdos  (Figure  6),  so  we  anticipated 
significant  correlations  with  parameteiizations  of  sediment  thickness.  However,  the  data 
display  relatively  weak  correlations  (absolute  values  <  0.38)  with  path-averaged  sediment 
thickness  (Figures  A2p,  A3p,  A4p,  A5p,  A6p),  and  only  slighUy  stronger  coirelations 
(absolute  values  <  0.41)  with  maximum  sediment  thickness  on  each  path  (Figures  A2r, 
A3r,  A4r,  A5r,  A6r).  Ratios  involving  Lg  tend  to  have  the  highest  correlations,  and  those 
are  in  the  lower  fiequency  ranges.  The  localized  nature  of  the  sedimentary  structures  in  the 
region,  which  thus  affect  only  subsets  of  raypaths  may  preclude  any  general  trends  from 
emerging  on  a  regional  basis,  and  the  low  correlations  do  not  indicate  that  sedimentary 

structures  are  not  important  on  a  case  by  case  basis. 

These  linear  regressions  suggest  that  important  waveguide  parameters  include  mean 
crustal  thickness,  mean  elevation  and/or  rms  roughness.  However,  these  parametric 
measures  are  not  independent,  and  in  fact,  they  all  exhibit  significant  correlations  with 
distance  that  must  be  allowed  for  before  any  physical  interpretation  or  optimal  path 
correction  strategy  can  be  advanced.  This  is  also  true  when  combinations  of  parameters  are 
made,  such  as  using  the  product  path  length  times  average  elevation,  which  gives 
particularly  high  correlations. 

Figure  9a  illustrates  the  coUinearity  of  some  of  the  waveguide  properties  with 
propagation  distance  for  our  data  set.  Particularly  strong  linear  correlations  are  found 
between  mean  elevation  and  distance  and  between  surface  roughness  and  distance.  Given 
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Figure  9.^  a)  Path  specific  values  of  mean  elevation,  nns  topographic  roughness,  mean 
thickness,  and  mean  sediment  thickness  as  functions  of  propagation  distance  from 
WMQ.  CC,  linear  coirelation  coefficient;  SD,  standard  deviation  of  the  linear  regression; 
K,  slope  of  Ae  regression  line.  The  dotted  lines  indicate  ±  1  standard  deviation  about  the 
regression  line,  b)  Correlations  between  mean  elevation  and  ims  topographic  roughness 
with  mean  crustal  Aickness  (top)  and  rms  topographic  roughness  and  rms  surface  slope 
with  mean  paA  elevation  (bottom).  CC,  linear  correlation  coefficient;  SD,  standard 
deviation  of  Ae  ^ear  regression;  K,  slope  of  Ae  regression  line.  The  dotted  lines  inAcate 
±  1  standard  deviation  about  Ae  regression  line. 

9.4 


MEAN  ELEVATION  (M/KM)  ROUGHNESS  (KM) 

-0.5  0.5  1.5  2.5  3.5  4.5  -0.5  0.5  1.5 


m 


MEAN  ELEVATION  (KM) 


O 


MEAN  ELEVATION  (KM) 


Figure  9.  a)  Path  specific  values  of  mean  elevation,  rms  topographic  roughness,  mean 
crustal  thickness,  and  mean  sediment  thickness  as  functions  of  propagation  distance  from 
WMQ.  CC,  linear  correlation  coefficient;  SD,  standard  deviation  of  the  linear  regression; 
K,  slope  of  the  regression  line.  The  dotted  lines  indicate  ±  1  standard  deviation  about  the 
regression  line,  b)  Correlations  between  mean  elevation  and  rms  topographic  roughness 
with  mean  crustal  thickness  (top)  and  rms  toi»graphic  roughness  and  rms  surface  slope 
with  mean  path  elevation  (bottom).  CC,  linear  correlation  coefficient;  SD,  smdard 
deviation  of  the  linear  regression;  K,  slope  of  the  regression  line.  The  dotted  lines  indicate 
±  1  standard  deviation  alwut  the  regression  line. 
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the  common  procedure  of  correcting  discriminants  for  distance,  the  issue  that  arises  is 
whether  or  not  there  is  any  value  in  correcting  for  more  than  distance.  Significant 
correlations  also  exist  between  various  parameters  of  waveguide  structure.  For  example, 
correlations  between  mean  crustal  thickness  and  mean  elevation  or  rms  roughness  are 
particularly  strong,  as  shown  in  Figure  9b.  This  reflects  isostatic  compensation  of  the 
crust,  as  do  surface  roughness  and  slope  correlations  with  mean  elevation  (Figure  9b). 
Previous  studies  have  detected  similar  correlations  between  parameters  (e.g.  Zhang  and 
Lay,  1994;  Zhang  et  al.,  1994),  which  makes  it  difficult  to  isolate  the  physical  effects 
controlling  regional  phase  energy  partitioning. 

However,  it  is  straightforward  to  test  whether  or  not  there  is  correlation  with  a 
given  parameter  after  correcting  for  the  effects  of  other  parameters.  We  did  this  for  all 
waveguide  parameters,  after  correcting  the  various  frequency  dependent  logarithmic 
amplitude  ratios  for  the  best  linear  dependence  on  distance.  The  “corrected”  data  thus 
correspond  to  those  values  that  are  typically  entered  into  a  discriminant  procedure.  If  any 
of  the  variance  in  the  residuals  can  be  reduced  by  correction  for  additional  parameters  we 
would  expect  to  see  the  discriminant  performance  improve.  We  find  that  the  process  of 
correcting  for  distance  systematically  reduces  the  correlations  with  other  parameters  (except 
in  a  few  cases  where  very  litde  correlation  initially  existed,  such  as  for  Pg/Lg  dependence 
on  mean  sediment  thickness),  typically  lowering  the  correlation  coefficients  by  about  0. 1  to 
0.3.  Figure  10  shows  the  regression  for  distance-corrected  Log(Pg/Lg)  measurements 

with  mean  crustal  thickness.  While  the  correlations  are  reduced  in  aU  frequency  bands 
relative  to  the  raw  data  correlations  in  Figure  8,  significant  correlation  remains  in  the  two 
lower  frequency  bands,  and  this  can  be  exploited  to  further  reduce  the  scatter  in  the  data  in 
an  objective  fashion.  The  Appendix  includes  results  for  distance  correction  of  aU  five 
ratios,  for  mean  elevation  (Figures  A2c,  A3c,  A4c,  A5c,  A6c),  roughness  (Figures  A2e, 
A3e,  A4e,  A5e,  A6e),  skewness  (Figures  A2i,  A3i,  A4i,  A5i,  A6i),  slope  (Figures  A2k, 
A3k,  A4k,  A5k,  A6k),  crustal  thickness  (Figures  A2m,  A3m,  A4m,  A5m,  A6m),  and 
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Figure  10.  Log(Pg/Lg)  as  a  function  of  mean  crustal  thickness  after  distance  correction  in 
various  frequency  bands.  CC,  linear  correlation  coefficient;  SD,  standard  deviation  of  the 
linear  regression;  K,  slope  of  the  regression  line.  The  dotted  lines  indicate  ±  1  standard 
deviation  about  the  regression  line. 
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mean  sediment  thickness  (Figures  A2q,  A3q,  A4q,  A5q,  A6q).  While  covariance  among 
the  path  properties  still  precludes  identifying  the  controlling  parameters,  the  two  low 
frequency  bands  indicate  empirical  dependence  (absolute  value  of  linear  correlations  >0.3, 
with  most  >  0.4  and  up  to  0.67)  on  mean  elevation,  surface  roughness,  and  crustal 
thickness,  and  weaker  dependence  (absolute  value  of  linear  correlations  >  0.2)  on 
skewness  (except  for  Pg/Sn  and  Sn/Lg),  and  mean  sediment  thickness  (only  for  Pg/Lg  and 
Sn/Lg). 

These  observations  make  a  solid  case  for  the  potential  value  of  empirical  correction 
for  path-specific  properties  beyond  simple  path  length.  In  order  to  develop  a  model  which 
explicitly  accounts  for  (or  at  least  allows  for)  collinearity  of  the  various  path  parameters,  we 
use  multivariate  regression  analysis.  Keeping  in  mind  that  our  objective  is  to  obtain  a 
stable  procedure  for  reducing  variance  in  the  observations,  we  do  not  need  to  uniquely 
isolate  the  causal  effects  between  coUinear  parameters,  as  long  as  their  trade-off  is 
incorporated  into  the  model. 


MULTIVARIATE  REGRESSION 


Multivariate  regression  allows  the  determination  of  the  best  set  of  predictor 
variables  from  a  set  of  parametric  path  properties  giving  an  empirical  model  minimizing  the 
variance  in  the  data.  An  initial  list  of  variables  is  defined,  based  on  the  preceding  separate 
correlation  analysis,  then  an  optimal  combination  is  determined.  Variables  with  such 
strong  collinearity  with  the  other  parameters  that  they  fail  to  contribute  significantly  to  the 
variance  reduction  are  excluded  from  the  ensuing  statistical  analysis  to  minimize  the  model 
complexity.  For  example,  maximum  and  minimum  crustal  thickness  and  maximum 
sediment  thickness  are  not  used  in  the  multivariate  regression  analysis  because  they  provide 
very  redundant  information  about  the  shape  of  the  Moho  discontinuity  and  basement 
structure  to  simple  (and  more  reliable)  path-averaged  values.  Surface  topographic  slope  is 
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excluded  because  it  almost  never  had  correlations  with  the  data.  We  retain  five  parameters 
to  characterize  the  path  properties:  x(l),  path  length;  x(2),  mean  altitude;  x(3),  rms 
roughness;  x(4),  mean  crustal  thickness;  and  x(5),  mean  sediment  thickness.  We  construct 
a  model  in  which  these  five  parameters  are  initially  assumed  to  influence  a  given  frequency 
dependent  regional  phase  ratio,  Y,  such  as  Log(Pg/Lg),  in  the  form 

Yj  =^0  A*U '*'^2^i2  ^  ~  2,  ...,  n,  (1) 

where  Y^  is  the  i-th  observation  that  is  assumed  to  be  a  linear  combination  of  the 
Xii,Xi2,...,Xi5  predictor  variables  with  linear  coefficients  ef  is  the  random 

error  associated  with  measurement  of  Yj,  and  n  is  the  number  of  observed  values.  We  do 
not  always  consider  all  5  parameters,  and  if  k  is  die  number  of  independent  variables 
considered,  the  model  has  m  =  k+1  unknown  parameters  (m  varies  up  to  6).  The  estimates 
of  pj  that  we  obtain  are  indicated  by  bj. 

It  is  desirable  to  follow  an  objective  procedure  for  limiting  the  number  of 
independent  parameters  while  still  achieving  near  optimal  variance  reduction  in  the  fit  to  the 
data.  A  variable  selection  technique  commonly  in  use  is  the  stepwise  regression  procedure, 
which  has  two  possible  versions;  forward  selection  and  backward  elimination.  In  our 
regression  model,  the  path  parameters  certainly  have  multicoUinearity,  which  is  mainly  a 
problem  for  uniquely  isolating  the  importance  of  each  parameter.  However,  this  is  not  all 
that  critical  in  this  case,  for  we  are  not  attempting  to  extract  a  deterministic  model  for  the 
medium.  For  example,  when  normal  distance  corrections  are  determined  there  is  no 
attempt  to  find  a  model  that  predicts  them,  and  the  same  coUinearity  with  other  path 
parameters  ensures  that  this  would  be  nonunique.  Empirical  path  corrections  need  not  be 
based  on  absolute  models,  and  hence  we  can  proceed  to  find  a  purely  empirical  model 
which  gives  favorable  variance  reduction  even  if  some  parameters  are  redundant  or  trade¬ 
off  in  their  contributions  to  the  fit 
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To  avoid  some  drawbacks  that  result  from  multicollinearity  when  applying  the 
stepwise  regression  technique,  we  did  not  use  backward  elimination  as  did  Zhang  et  al. 
(1994).  Instead,  we  compute  all  possible  linear  regressions  for  our  combination  of  5  path 
parameters  and  then  select  a  preferred  set  of  best  parameters  for  each  data  type  based  on 
statistical  properties  of  the  solutions.  We  use  two  common  criteria:  (1)  the  mean  square 
error  (MSE)  criterion,  and  (2)  the  Cp-statistic  criterion.  The  MSE  is  the  residual  variance: 


MSE 


SSE 


£(yi-yi)^ 


n-m  n-m 


(2) 


where  SSE  represents  the  sum  of  squared  residuals,  y,  are  observed  values,  and  yj  are 
predicted  values  for  the  model.  The  MSE  criterion  minimizes  the  variance,  taking  into 
account  the  number  of  parameters  in  the  model  (m),  but  it  does  not  assess  the  resolution  or 
contribution  of  each  parameter.  Mallow’s  Cp-statistic  criterion  is  given  by: 

Cp=^-(n-2p)  (3) 

where  SSEp  is  the  sum  of  squared  errors  for  the  model  obtained  for  a  subset  of  p  variables 

out  of  a  total  of  k  predictor  variables  being  considered;  and  S^  =  MSE(Xj,X2,...,x^)  is  the 
estimation  variance  based  on  all  k  predictor  variables.  After  aU  possible  regression 
combinations  are  determined,  Cp  is  computed  for  each  case,  and  regressions  that  have  Cp 

values  close  to  p  (i.e.  low  values)  are  the  most  desirable. 

Table  1  lists  the  model  coefficients  found  for  all  possible  regressions  for 
Log(Pg/Lg)  data  in  the  frequency  band  of  0.75  to  1.5  Hz.  Table  2  gives  goodness  of  fit 
criteria  for  all  possible  regressions  for  the  same  data  set  Table  2  indicates  that  there  are  two 
models  that  stand  out;  the  models  with  parameters  (Xi,X2,X4,X5)  or  (Xi,X2,X3,X4,X5). 
Both  models  have  high  multiple  coefficient  of  determination  (R^),  small  residual  variances, 
and  small  Cp  values.  Further  analysis  indicates  that  the  second  model  containing  all  five 
waveguide  parameters  is  not  particularly  appealing  because  the  estimate  of  the  regression 
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TABLE  1.  Estimated  regression  coefficients  for  0.75-1.5  Hz  Log(Pg/Lg)  data  observed  at  WMQ 


Predictor 

variables 

bo 

bi 

b2 

b3 

b4 

b5 

-0.114 

xi  (Distance) 

-0.881 

0.000596 

x2(mH) 

-1.218 

0.476 

x3  (Roughness) 

-1.161 

0.895 

x4  (mMoho) 

-5.458 

0.101 

x5  (mSed.) 

-0.086 

-0.0123 

XI,  x2 

-1.339 

0.000254 

0.387 

XI,  x3 

-1.345 

0.000281 

0.744 

XI,  x4 

-5.090 

0.000367 

0.0851 

XI,  x5 

-1.174 

0.000692 

0.0748 

x2,x3 

-1.425 

0.267 

0.592 

x2,x4 

-4.234 

0.257 

0.666 

x2,x5 

-1.500 

0.524 

0.0742 

X3,X4 

-3.796 

0.510 

0.0583 

X3,x5 

-1.143 

0.894 

-0.00754 

x4,  x5 

-5.607 

0.103 

0.0268 

XI,  x2,  x3 

-1.490 

0.000165 

0.223 

0.553 

XI,  x2,  x4 

-4.438 

0.000274 

0.157 

0.0682 

XI,  x2,  X5 

-1.774 

0.000356 

0.419 

0.102 

XI,  x3,  x4 

-4.038 

0.000290 

0.346 

0.0595 

XI,  x3,  x5 

-1.461 

0.000336 

0.717 

0.0338 

XI,  x4,  x5 

-5.430 

0.000467 

0.0857 

0.0791 

x2,  x3,  x4 

-3.209 

0.202 

0.395 

0.0409 

X2,  x3,  x5 

-1.567 

0.311 

0.544 

0.0419 

x2,  x4,  x5 

-4.341 

0.0308 

0.0637 

0.0628 

x3,  x4,  x5 

-3.902 

0.500 

0.0600 

0.0132 

XI,  x2,  X3,  x4 

-3.590 

0.000218 

0.134 

0.311 

0.0477 

XI,  x2,  X3,  x5 

-1.746 

0.000246 

0.271 

0.459 

0.0659 

XI,  x2,  x4,  x5 

-4.658 

0.000346 

0.197 

0.0645 

0.0907 

XI,  x3,  x4,  x5 

-4.646 

0.000395 

0.237 

0.0680 

0.0646 

x2,  x3,  x4,  X5 

-3.451 

0.249 

0.332 

0.0428 

0.0468 

XI,  x2,  x3,  x4. 

-4.185 

0.000324 

0.180 

0.164 

0.0541 

0.0796 

x5 


TABLE  2.  Goodness  of  fit  criteria  for  regressions  for  0.75-1.5  Hz  Log(Pg/Lg)  data  WMQ 


Predictor  Variables 

r2 

SSE 

VARIANCE 

Cp 

- 

19.129 

0.226 

XI  (Distance) 

0.314 

13.114 

0.154 

147.070 

x2(mH) 

0.508 

9.414 

0.111 

82.158 

x3  (Roughness) 

0.559 

8.444 

0.099 

65.140 

x4  (iiiM(±o) 

0.561 

8.393 

0.099 

64.246 

X5  (mSed.) 

0.001 

19.103 

0.225 

252.140 

XI,  X2 

0.547 

8.658 

0.103 

70.895 

XI,  x3 

0.613 

7.412 

0.088 

49.035 

XI,  x4 

0.667 

6.376 

0.076 

30.860 

XI,  x5 

0.357 

12.309 

0.147 

134.947 

x2,  X3 

0.654 

6.620 

0.079 

35.140 

x2,  x4 

0.645 

6.794 

0.081 

38.193 

X2,  x5 

0.552 

8.573 

0.102 

69.404 

X3,  x4 

0.642 

6.847 

0.082 

39.123 

X3,X5 

0.559 

8.435 

0.100 

66.982 

x4,  X5 

0.568 

8.273 

0.098 

64.140 

XI,  X2,  x3 

0.670 

6.313 

0.076 

31.754 

XI,  X2,  x4 

0.691 

5.914 

0.071 

24.754 

XI,  X2,  X5 

0.623 

7.205 

0.087 

47.404 

XI,  x3,  x4 

0.699 

5.753 

0.069 

21.930 

XI,  x3,  X5 

0.621 

7.256 

0.087 

48.298 

XI,  x4,  x5 

0.714 

5.477 

0.066 

17.088 

X2,  X3,  X4 

0.689 

5.940 

0.072 

25.211 

x2,  x3,  X5 

0.667 

6.370 

0.077 

32.754 

X2,  x4,  x5 

0.676 

6.197 

0.075 

29.719 

X3,  x4,  x5 

0.644 

6.819 

0.082 

40.632 

XI,  X2,  x3,  x4 

0.717 

5.422 

0.066 

18.123 

XI,  X2,  x3,  x5 

0.698 

5.769 

0.070 

24.211 

XI,  x2,  x4,  x5 

0.751 

4.763 

0.058 

6.579 

XI,  x3,  X4,  x5 

0.727 

5.214 

0.064 

14.474 

x2,  x3,  X4,  x5 

0.706 

5.630 

0.069 

21.772 

XI,  x2,  X3,  x4,  X5 

0.757 

4.644 

0.057 

6.474 
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coefficient  for  the  x(3)  term  has  a  standard  deviation  which  is  quite  large  (0.164  ±  0.1 13), 
and  the  partial  F-statistic  for  this  coefficient  F(jSj  =  0)  value  is  2.10,  smalls  than  2.77,  the 

10%  light-tail  point  of  F(0.90,  1,  81),  indicating  that  we  cannot  reject  at  the  90% 
confidence  level  the  null  hypothesis  that  pj  =  0.  In  contrast,  the  four-parameter  model  has 
satisfactory  accuracy  in  the  estimates  of  all  model  coefficients  and  passed  an  F-test  of  tiie 
overall  regression  model,  for  whether  at  least  one  of  the  coefficients  is  zero.  Table  3  lists 
the  model  parameters  for  this  case.  All  partial  F(Pj  =  0)  values  are  greater  than  6.96,  the 
1%  right-tail  point  of  F(0.99,  1,  82),  thus  the  model  passed  all  partial  F-tests,  indicating 
that  the  null  hypothesis  that  pj  =  0  can  be  rejected  with  99%  confidence.  In  addition,  aU  t- 
statistics  are  substantially  greater  than  rejection  points  at  99.9%  confidence  level,  indicating 
the  Pj  is  nonzero  and  the  four  model  parameters  are  related  to  the  amplitude  ratio  being 
considered.  If  we  perform  a  formal  stepwise  backward  elimination  we  obtain  the  same 
model.  Thus,  a  model  which  includes  the  four  path  parametCTS,  distance,  mean  altitude, 
mean  crustal  thickness,  and  mean  sediment  thickness,  is  preferred.  This  model  explains 
75.1%  of  the  total  variation  in  the  87  amplitude  Log(Pg/Lg)  measurements  for  the 

frequency  band  0.75  to  1.5  Hz. 

The  primary  objective  here  is  to  reduce  the  scatter  in  the  data,  and  it  is  helpful  to 
visualize  the  reduction  that  occurs  for  various  parameter  combinations.  Figure  1 1  displays 
the  scatter  in  the  raw  and  corrected  data  for  various  models  obtained  by  multiple  regression. 
The  observations  are  plotted  as  a  function  of  magnitude,  and  ideally  there  should  be  no 
magnitude  dependence  (none  is  observed)  and  as  little  range  of  observations  as  possible  (so 
that  an  outlying  explosion  signal  could  be  detected).  Noting  the  large  scatter  in  the  raw 
measurements,  it  is  clear  that  for  this  low  fiequency  band,  effective  source  discrimination 
will  be  very  problematic  without  path  corrections.  Standard  distance  correction  gives  a 
significant,  31%  variance  reduction,  but  this  is  actually  smaller  than  that  achieved  by 
correcting  for  just  mean  path  elevation  (50%  variance  reduction),  or  mean  crustal  thickness 
(65%  reduction).  The  added  value  of  multiple  parameter  regressions  is  also  apparent,  with 
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Table  3.  The  best  set  of  predictor  variables  for  0.75-1.5  Hz  Pg/Lg  data  at  WMQ 


Variable 

Estimated  bj 

Standard  Deviation 

F(Pj=0) 

t(Pi=0) 

x(l) 

0.000364 

0.0000733 

24.7 

4.97 

x(2) 

0.197 

0.0562 

12.3 

3.51 

x(4) 

0.0645 

0.00996 

42.0 

6.48 

x(5) 

Intercept 

0.0907 

-4.658 

0.0204 

19.8 

4.45 

r2=  0.751  =  0.739 

The  Total  Sum  of  Square  (SST)  =  19.129 

The  Error  Sum  of  Square  (SSE)  =  4.763 

The  Error  Mean  Square  (MSE)  =  0.058 

Standard  Deviation  =  0.241 

Number  of  Observation=  87  F(  4,  82)  =  61.834 
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Figure  11.  Log(Pg/Lg)  measurements  in  the  frequency  band  of  0.75  to  1.5  Hz  as  a 
function  of  event  magnitude  with  various  level  of  propagation  correction.  Orig.,  observed 
data;  Dist.,  corrected  for  path  length;  mH,  corrected  for  mean  altitude  on  each  path;  mM, 
corrected  for  mean  crustal  thickness  on  each  path;  Dist.+Rough.,  corrected  for  a  two- 
parameter  model  involving  distance  and  ims  roughness;  Dist.+mM,  corrected  for  a  two- 
parameter  model  involving  distance  and  mean  crustal  thickness;  Dist.+mM+mS,  corrected 
for  a  three  parameter  model  involving  distance,  mean  crustal  thickness,  and  mean  sediment 
thickness;  Dist.+mH+mM+mS,  corrected  for  a  four  parameter  model  involving  distance, 
mean  altitude,  mean  crustal  thickness,  and  mean  sediment  thickness. .  In  each  case  the 
variance  of  the  observations  is  indicated  below  the  plot  by  Var. 
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Corrected  LOG[Pg/Lg]  Corrected  LOG[Pg/Lg] 


two  parameter  models  doubling  the  variance  reduction  relative  to  simple  distance 
corrections,  while  the  preferred  four  parameter  model  gives  74%  variance  reduction,  more 
than  twice  that  of  the  standard  distance  correction.  Figure  12  portrays  the  same  results  as 
histograms  of  observations,  corrected  for  various  models.  The  initial  data  distribution 
evolves  toward  an  increasingly  symmetric,  normal  distribution  with  decreasing  standard 
deviation  as  the  model  complexity  increases. 

Qualitatively  similar  results  are  found  for  the  lowest  frequency  band  for  other 
phase  ratios,  as  well  as  for  frequencies  from  1.5  to  3.0  Hz.  The  results  of  multiple 
regressions  for  the  Pg/Lg  observations  in  the  1.5  to  3.0  Hz  band  are  summarized  in  Table 
4,  and  the  details  of  the  best  model  are  given  in  Table  5  (results  for  other  phase  pairs  will 
be  presented  in  a  separate  paper).  Similar  to  the  lower  frequency  band,  distance  correction 
alone  reduces  the  variance  by  42%,  while  mean  path  elevation  and  mean  roughness 
provide  better  variance  reduction  on  their  own.  A  4-parameter  model  involving  distance, 
mean  elevation,  surface  roughness  and  mean  sedimentary  layer  thickness  provides  the  best 
overall  variance  reduction  (70%),  with  all  coefficients  passing  the  partial  F(Pj  =  0)  and  t- 
tests  (Table  5).  While  the  overall  variance  reduction  from  the  additional  path  parameters  is 
reduced  somewhat  in  tiiis  band  relative  to  the  lowest  frequency  band  (as  expected,  given 
the  decreasing  correlations  observed  as  frequency  increases  in  Figures  7,  8  and  10),  the 
improvement  relative  to  distance  correction  alone  is  still  significant.  The  residual  variance 
of  the  observations  for  this  model  is  74%  of  that  for  the  best  fitting  model  in  the  0.75- 1 .5 
Hz  band. 

Higher  frequency  (>3  Hz)  observations  exhibit  weaker  dependence  on  path 
properties  for  all  of  the  phase  ratios,  but  path  correction  procedures  can  still  reduce  the 
scatter,  improving  discrimination.  The  results  of  multiple  regressions  for  the  Pg/Lg 
observations  in  the  3.0-6.0  Hz  band  are  shown  in  Table  6.  In  this  case,  distance  is  the  best 
single  parameter  for  correction,  giving  a  20%  variance  reduction,  but  a  multiple-regression 
model  involving  distance,  mean  elevation,  mean  roughness,  and  mean  crustal 
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TABLE  4.  Goodness  of  fit  criteria  for  regressions  for  1. 5-3.0  Hz  Log(Pg/Lg)  data  at  WMQ. 


Predictor  Variables 

R2 

SSE 

Variance 

Cp 

- 

12.171 

0.1415 

XI  (Distance) 

0.431 

6.928 

0.0815 

77.742 

X2(mH) 

0.477 

6.360 

0.0748 

65.563 

X3  (Roughness) 

0.528 

5.747 

0.0676 

50.341 

X4  (mMoho) 

0.360 

7.789 

0.0916 

98.719 

X5  (mSed.) 

0.006 

12.094 

0.1423 

197.603 

XI,  X2 

0-585 

5.050 

0.0601 

36.169 

XI,  X3 

0.652 

4.236 

0.0504 

17.283 

XI,  X4 

0.591 

4.979 

0.0593 

34.522 

XI,  X5 

0.471 

6.438 

0.0766 

68.374 

X2,X3 

0.617 

4.667 

0.0556 

27.283 

X2,X4 

0.514 

5.919 

0.0705 

56.332 

X2,X5 

0.502 

6.059 

0.0721 

59.580 

X3,X4 

0.536 

5.651 

0.0673 

50.114 

X3,  X5 

0.532 

5.695 

0.0678 

51.135 

X4,  X5 

0.360 

7.787 

0.0927 

99.673 

XI,  X2,  X3 

0.684 

3.842 

0.0463 

10.142 

XI,  X2,  X4 

0.627 

4.542 

0.0547 

26.383 

XI,  X2,  X5 

0.652 

4.237 

0.0510 

19.306 

XI,  X3,  X4 

0.661 

4.121 

0.0496 

16.615 

XI,  X3,  X5 

0.662 

4.109 

0.0495 

16.336 

XI,  M,  X5 

0.634 

4.451 

0.0536 

24.271 

X2.  X3,  X4 

0.617 

4.661 

0.0562 

29.143 

X2,  X3,  X5 

0.620 

4.620 

0.0557 

28.193 

X2,  X4,  X5 

0.533 

5.679 

0.0684 

52.763 

X3,  X4,  X5 

0.538 

5.623 

0.0677 

51.464 

XI,  X2,  X3,  X4 

0.685 

3.833 

0.0467 

11.933 

XI,  X2,  X3,  X5 

0.711 

3.520 

0.0429 

4.671 

XI,  X2,  X4,  X5 

0.685 

3.838 

0.0467 

12.049 

XI,  X3,  X4,  X5 

0.678 

3.915 

0.0477 

13.835 

X2,  X3,  X4,  X5 

0.621 

4.616 

0.0563 

30.100 

XI,  X2,  X3,  X4,  X5 

0.713 

3.489 

0.0431 

5.951 
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Table  5.  The  best  set  of  predictor  variables  for  1.5-3.0  Hz  Log(Pg/Lg)  data  at  WMQ. 


Variable 

bj 

Standard  Deviation 

F(bi=0) 

t(bi=0) 

x(l) 

0.000333 

0.0000658 

25.6 

5.06 

x(2) 

0.171 

0.0461 

13.7 

3.70 

x(3) 

0.324 

0.0794 

16.7 

4.09 

x(5) 

0.0507 

0.0185 

7.49 

2.74 

Litercept 

-1.226 

r2=  0.711  =  0.697 

The  Total  Sum  of  Square  (SST)=  12.171 

The  Error  Sum  of  Square  (SSE)=  3.520 

The  Error  Mean  Square  (MSE)=  0.043 

Standard  Deviation=  0.207 

Number  of  Observation=  87  F(  4, 82)=  50.375 
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TABLE  6.  Goodness  of  fitcriteria  for  regressions  for  3.0-6.0  Hz  Log(Pg/Lg)  data  at  WMQ. 


Predictor  Variables 

R2 

SSE 

Variance 

Cp 

- 

5.489 

0.0638 

XI  (Distance) 

0.207 

4.354 

0.0512 

20.667 

X2(mH) 

0.162 

4.601 

0.0541 

26.548 

X3  (Roughness) 

0.168 

4.568 

0.0537 

25.762 

X4  (mMoho) 

0.012 

5.420 

0.0638 

46.048 

X5  (mSed.) 

0.011 

5.430 

0.0639 

46.286 

XI,  X2 

0.239 

4.176 

0.0497 

18.428 

XI,  X3 

0.255 

4.090 

0.0487 

16.381 

XI,  X4 

0.209 

4.342 

0.0517 

22.381 

XI,  X5 

0.214 

4.313 

0.0513 

21.690 

X2.X3 

0.202 

4.382 

0.0522 

23.333 

X2,X4 

0.204 

4.367 

0.0520 

22.976 

X2,  X5 

0.163 

4.596 

0.0547 

28.428 

X3.X4 

0.252 

4.108 

0.0489 

16.810 

X3.X5 

0.177 

4.517 

0.0538 

26.548 

X4.X5 

0.020 

5.378 

0.0640 

47.048 

XI,  X2,  X3 

0.261 

4.055 

0.0489 

17.548 

XI,  X2,  X4 

0.277 

3.967 

0.0478 

15.452 

XI,  X2,  X5 

0.252 

4.107 

0.0495 

18.786 

XI,  X3,  X4 

0.335 

3.652 

0.0440 

7.952 

XI,  X3,  X5 

0.256 

4.082 

0.0492 

18.190 

XI,  X4,  X5 

0.216 

4.302 

0.0518 

23.429 

X2,  X3,  X4 

0.349 

3.571 

0.0430 

6.024 

X2,  X3,  X5 

0.202 

4.379 

0.0528 

25.262 

X2,  X4,  X5 

0.207 

4.354 

0.0525 

24.667 

X3,  X4,  X5 

0.277 

3.968 

0.0478 

15.476 

XI,  X2,  X3,  X4 

0.380 

3.402 

0.0415 

4.000 

XI,  X2,  X3,  X5 

0.266 

4.031 

0.0492 

18.976 

XI,  X2,  X4,  X5 

0.294 

3.876 

0.0473 

15.286 

XI,  X3,  X4,  X5 

0.335 

3.647 

0.0445 

9.833 

X2,  X3,  X4,  X5 

0.352 

3.559 

0.0434 

7.738 

xi.X2,X3,  X4,  X5 

0.380 

3.400 

0.0420 

5.952 
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thickness  provides  a  35%  variance  reduction,  with  aU  coefficients  being  statisticaUy 
dgnifirant  at  a  99%  confidence  level.  It  is  particularly  notable  that  the  residual  variance 
for  the  preferred  4-parameter  model  is  essentially  identical  to  that  for  the  preferred  4- 
parameter  model  for  the  1.5-3  Hz  band.  This  indicates  that  we  are  indeed  achieving  the 
optimal  reduction  of  variance  in  the  lower  firequency  band,  and  discrimination  may  be 
equally  good  across  the  passband  for  diis  data  set.  Comparable  results  are  found  for  the 
other  ratios  (to  be  presented  elsewhere).  The  6-9  Hz  band  displays  very  Utde  dependence 
on  path  parameters,  and  we  do  not  show  results  for  that  case  as  we  believe  that  our  signal- 
to-noise  ratios  are  marginal.  It  may  be,  however,  that  when  data  sets  for  more  localized 
source  regions  are  considered  both  distance  and  other  effects  may  be  significant  at  high 
frequency.  Alternatively,  in  some  regions,  no  corrections  may  be  warranted,  depending 
on  the  nature  of  the  regional  propagation.  Our  recommendation  is  that  every  station/region 
be  treated  on  a  case  by  case  basis,  with  multiple  regression  analysis  be  explored  in  each 
case  and  the  resulting  used  when  the  model  reduces  the  variance  favorably.  Modeling 
efforts  need  to  be  undertaken  to  place  die  empirical  correlations  on  a  sound  physical  basis. 

DISCUSSION  AND  CONCLUSIONS 

While  the  reduction  in  variance  associated  with  the  multiple  regression  analysis  for 
Pg/Lg  phases  recorded  at  WMQ  is  very  encouraging,  it  must  be  emphasized  that  die 
corresponding  model  is  empirical,  is  unlikely  to  generalize  because  of  die  extreme  path 
complexity  in  the  region,  and  is  only  useful  for  the  fiequency  bands  lower  than  3  Hz.  The 
empirical  nature  of  the  model  is  not  very  troubling,  as  the  same  can  be  said  for  all  standard 
distance-correction  procedures,  although  they  are  often  motivated  by  quantitative  wavefield 
behavior  in  one-dimensional  models.  The  other  path  attributes  considered  explicidy  require 
three-dimensional  modeling  which  has  not  been  performed  extensively,  but  the  general 
tendencies  are  consistent  with  past  empirical  and  two-dimensional  modeling  efforts.  The 
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fact  that  the  path  complexity  around  WMQ  is  very  acute,  particularly  for  partial  and  total 
blockage  of  Lg  phases  from  events  in  the  Tibetan  Plateau,  is  a  realistic  fact  of  life  for  CTBT 
monitoring.  We  often  will  have  only  a  few  stations  in  tectonically  complex  regions  and 
must  have  procedures  that  allow  for  great  path  complexity.  However,  we  do  expect  that 
different  waveguide  parameters  may  emerge  as  most  important  in  regions  of  less  crustal 
heterogeneity,  and  in  some  cases  standard  correction  for  path  length  may  be  all  that  is 
warranted.  The  lack  of  strong  path  (or  distance)  dependence  of  die  higher  frequency  ratios 
is  interesting,  as  it  has  commonly  been  found  that  discrimination  is  poor  at  low  frequencies 
and  improves  above  3  Hz.  Improving  the  performance  at  low  frequencies  is  thus  highly 
desirable,  and  factor  of  2  to  3  reductions  in  the  scatter  of  discriminant  measures  for  an 
earthquake  population  should  help  directly.  Again,  this  is  particularly  important  for 
tectonically  active  areas,  whae  strong  attenuation  may  constrain  the  observed  bandwidth  of 
regional  phases  to  lower  frequencies. 

All  of  the  regional  phase  amplitude  ratios  considered;  Pn/Lg,  Pg/Lg,  Pn/Sn»  Pg/Sn» 
and  Sn/Lg,  exhibit  frequency  dependent  variations  with  propagation  distance,  and 
correction  for  this  effect  is  necessary  for  discrimination.  However,  as  we  have  shown,  the 
uncorrected  data  sometimes  correlate  better  with  specific  path  waveguide  properties  than 
they  do  with  distance.  This  leaves  open  the  issue  of  what  are  the  controlling  effects  on 
regional  phases,  with  factors  that  are  probably  important  being  geometric  spreading, 
attenuation,  scattering,  sedimentary  basins,  changes  in  Moho  depth,  rough  surface 
topography,  and  internal  crustal  structure.  There  is  a  need  for  extensive  modeling  of 
realistic  (3-dimensional)  crustal  waveguides  with  structure  at  a  wide  variety  of  scale 
lengths,  if  we  are  to  fully  understand  the  nature  of  these  regional  phases. 

The  WMQ  data  exhibit  significant  Lg  blockage,  and  much  of  the  value  of  correction 

for  path-specific  wave  guide  properties  stems  from  effects  of  the  Tibetan  Plateau.  Current 
discrimmant  calibration  procedures  are  not  well  defined,  but  appear  to  either  screen  data 
based  on  defined  regions  of  phase  blockage  (hoping  that  some  other  discriminant  can  be 
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applied  in  those  regions),  or  to  develop  separate  distance  relations  in  regionalized  areas  that 
have  baseline  shifts  in  their  regional  phase  energy  partitioning.  The  systematic 
relationships  observed  in  this  study  between  distance-corrected  regional  discrintinants  and 
parameters  of  the  wave  guide  structure  suggests  that  blockage  effects  involve  a  continuum 
of  effects  ratiier  than  bimodal  behavior.  A  regional  discriminant  calibration  strategy  that 
invokes  multivariate  regression  analysis  of  several  path  properties  such  as  mean  crustal 
thickness  and  mean  altitude  allows  all  data  at  a  given  station  in  a  complex  region  to  be 
calibrated  uniformly,  without  appealing  to  separate,  possibly  poorly  constrained,  distance 
corrections  for  discrete  subregions.  This  can  readily  be  implemented  in  a  CTBT 
monitoring  system.  The  extent  to  which  it  provides  high  confidence  discrimination  needs 
to  be  further  tested,  but  the  basic  strategy  appears  promising. 
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Abstract 


Regional  phase  amplitude  ratio  measurements  are  useful  for  discriminating  earth¬ 
quakes  from  explosions  and  will  be  essential  for  seismic  monitoring  of  a  Comprehen 
sive  Test  Ban  Treaty  (CTBT).  However,  such  measurements  have  been  made  in  vari¬ 
ous  ways,  complicating  comparison  of  different  studies.  We  evaluate  several  regional 
phase  amplitude  ratio  measurement  procedures.  Amplitude  measurements  are  made  on 
both  original  velocity-proportional  broadband  seismograms  and  instrument-deconvolved 
displacement  seismograms  in  the  time  and  frequency  domains.  Pn/Lg  and  Pg/Lg  ratios 
are  measured  on  vertical  component  waveforms  for  regional  earthquakes  observed  at 
the  Global  Seismic  Network  (GSN)  station  ABKT  (Alibek,  Turkmenistan).  Time 
domain  amplitude  measurements  are  made  for  narrow  band  filtered  waveforms  using 
signal  energy  measures  given  by  the  absolute  mean,  rms,  and  envelope  mean,  and  peak 
measures  given  by  the  absolute  maximum  and  envelope  maximum.  Frequency  domain 
measurements  are  made  by  computing  the  log-mean  amplitude  within  a  narrow  band 
from  the  broadband  spectrum  of  each  windowed  phase. 

Time  domain  amplitude  ratios  for  raw  velocity  and  instrument  deconvolved  dis¬ 
placement  seismograms  are  very  similar  (linear  correlations  ^  0.97),  indicating  that  the 
instrument  response  does  not  affect  the  time  domain  measurements  for  the  broadband 
GSN  data.  Time  domain  amplitude  ratios  made  using  the  energy  or  peak  measures  are 
well  correlated  among  themselves  (linear  correlations  >  0.97).  However,  the  energy 
and  peak  measures  are  slightly  less  well  correlated  with  each  other  (linear  correlations 


between  0.87-0.99).  For  the  time  domain  measures,  the  correlations  generally  degrade 
and  the  scatter  increases  for  the  higher  frequency  bands.  Time  and  frequency  domain 
measurements  are  well  correlated  for  the  lowest  frequency  band  (0.75-1.5  Hz).  How¬ 
ever,  as  the  frequency  band  increases,  the  correlation  decreases  and  the  slope  deviates 
from  1.0.  Time  domain  amplitude  ratios  (for  both  Pn/Lg  and  Pg/Lg)  are  consistently 
larger  than  the  frequency  domain  amplitude  ratios.  Investigation  of  the  frequency 
spectra  of  broadband  and  narrow  band  filtered  phases  revealed  that  significant  energy 
from  outside  the  pass  band,  particularly  from  lower  frequencies,  possibly  biases  the 
time  domain  measurements.  For  the  highest  frequency  pass  band  (6.0-9.0  Hz),  decon¬ 
volution  of  the  instrument  response  causes  differences  between  the  broadband  and 
filtered  spectra.  Better  agreement  between  time  and  frequency  domain  measurements 
is  found  when  linear-averaging  of  the  spectrum  is  used.  These  observations  suggest 
that  log  averaged  frequency  domain  measurements  are  favorable  over  time  domain 
measurements  because:  1)  frequency  domain  measurements  offer  complete  control  of 
which  frequencies  are  used;  and  2)  log-averaging  of  the  spectra  does  not  over-weight 
spectral  amplitudes  associated  with  lower  frequencies  within  the  pass  band.  Further¬ 
more,  spectral  measurements  are  computationally  more  efficient  and  convenient  for 
automation  than  time  domain  measurements. 

Introduction 

Regional  broadband  seismic  data  have  proven  to  be  effective  for  discriminating 
explosions  from  natural  earthquakes.  Discrimination  at  regional  distance  is  especially 
critical  for  smaller  events  (mj,  <  4.5)  as  these  events  may  not  be  detected  and  discrim¬ 
inated  teleseismically  (Blandford,  1981;  Pomeroy  et  al.,  1982).  Amplitude  and  spectral 
ratios  of  regional  phases,  particularly  Pn,  Pg  and  Lg,  have  been  shown  to  separate 
explosion  and  earthquake  populations  in  many  cases  (e.g.  Taylor  et  al.,  1988;  Walter  et 
al.,  1995;  Hartse  et  al.,  1996).  Complicating  the  situation,  however,  are  path  effects. 
Regional  variations  in  geologic,  tectonic  and  topographic  structure  lead  to  great  varia¬ 
bility  in  regional  phase  amplitudes.  Furthermore,  recordings  of  regional  events  can  be 
noisy  due  to  poor  signal  excitation  and  attenuation,  especially  at  high  frequencies  and 
distances  greater  than  several  hundred  kilometers.  Nonetheless,  measurements  of 
regional  phase  amplitude  ratios  can  provide  a  valuable  tool  for  identifying  small 
seismic  events  under  the  Comprehensive  Test  Ban  Treaty  (CTBT)  provided  that  varia¬ 
tions  in  earthquake  generated  regional  phases  are  determined  for  each  monitoring  sta¬ 
tion  and  region.  Efforts  are  under  way  to  map  variations  in  regional  phase  amplitude 
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ratios  and  to  understand  the  path  effects  in  support  of  CTBT  monitoring. 

Complicating  comparison  of  various  results,  however,  is  the  fact  that  regional 
phase  amplitude  measurements  have  been  made  in  a  number  of  ways,  including  tech¬ 
niques  operating  in  either  the  time  or  frequency  domain.  Time  domain  measures  have 
been  made  on  filtered  and  windowed  phases  such  as  peak  amplitudes,  rms  amplitudes 
and  envelope  measures  (e.g.  Nuttli,  1980;  Zhang  and  Lay,  1994;  Hartse,  et  al.,  1996). 
Alternatively,  frequency  domain  measures  have  been  used  (e.g.  Taylor  et  al.,  1989; 
Walter  et  al.,  1995).  Also,  time  domain  measurements  can  be  made  on  raw  velocity  or 
instrument  deconvolved  displacement  seismograms.  Thus,  it  is  important  to  know  if 
results  obtained  from  different  measurement  methodologies  agree  or  not  and,  if  they  do 
not,  to  understand  why. 

In  this  study  we  investigated  the  process  of  making  regional  phase  amplitude  ratio 
measurements.  Pn/Lg  and  Pg/Lg  amplitude  ratios  were  computed  using  both  time  and 
frequency  domain  methods  on  raw  velocity  and  instrument  deconvolved  displacement 
seismograms.  Results  were  compared  and  the  cause  of  differences  in  the  measurements 
were  investigated. 

Regional  Earthquake  Data  at  ABKT 

The  data  used  in  this  study  were  recorded  at  the  Global  Seismic  Network  (GSN) 
station  ABKT  in  Alibek,  Turkmenistan.  This  station  has  a  Streckeisen  STS-2  broad¬ 
band  seismometer  with  flat  velocity  response  from  approximately  0.01-10  Hz  (sample 
rate  =  20/sec).  We  collected  waveform  data  recorded  at  ABKT  for  earthquakes  in  the 
region  reported  by  the  National  Earthquake  Information  Service-Preliminary  Determi¬ 
nation  of  Epicenters  (NEIS-PDE)  for  the  years  1994  and  1995.  Reported  depths  were 
less  than  50  km  and  body-wave  magnitudes,  m^,  were  greater  than  or  equal  to  4.0. 
The  events,  paths  and  region  studied  are  shown  in  Figure  1. 

All  waveform  data  were  previewed  and  a  Pn  arrival  was  picked.  Noisy  data  were 
discarded.  For  this  study  only  the  vertical  component  was  considered.  The  instrument 
response  was  removed  and  the  data  were  integrated  to  displacement.  Regional  phases 
were  isolated  with  group  velocity  windows:  Pn  8.0  to  6.2  km/sec;  Pg  6.2  to  5.0 
km/sec;  and  Lg  3.6  to  3.0  km/sec.  In  order  to  reduce  possible  biases  due  to  poor  event 
locations  and  timing  errors,  all  group  velocity  windows  were  shifted  such  that  Pn 
arrives  at  7.9  km/sec.  Noise  was  sampled  in  a  35  sec  window  ending  5  sec  before  the 
Pn  arrival.  Figure  2  shows  an  example  of  the  phase  windowing. 


ABKT.BHZ  1994  54  11  54  mb=  5.3  Distance=  823  km  Filtered  0.75-9.0  Hz 


Propagation  of  regional  phases  in  this  region  is  strongly  influenced  by  the  com¬ 
plex  topographic  and  tectonic  structure.  High  elevations  of  the  Turkish-Iranian  Plateau 
are  associated  with  recent  volcanic  activity  and  low  sub-Moho  P-wave  velocities 
(Hearn  and  Ni,  1994).  Regions  of  low  Pn  velocities  are  strongly  correlated  with  a  zone 
of  inefficient  Sn  propagation  (Kadinsky-Cade  et  al.,  1981;  Rodgers  et  al.,  1997).  Low 
Q  within  the  crust  of  the  Iranian  Plateau  leads  to  inefficient  Lg  propagation  (Nuttli, 
1980;  Wu  et  al.,  1996).  The  topographic  relief  associated  with  the  boundaries  of  the 
Turkish-Iranian  Plateau  has  been  shown  to  weaken  and/or  block  propagation  of  the  Lg 
phase  (Kadinsky-Cade  et  al.,  1981;  Rodgers  et  al.,  1997).  Figure  3  illustrates  the  vari¬ 
ability  in  Lg  propagation  for  the  region.  Pg/Lg  amplitude  ratios  were  made  as 
described  below  using  the  time  domain  rms  amplitudes  in  four  frequency  bands.  The 
Pg/Lg  ratio  for  each  path  is  plotted  and  color  coded.  Note  that  the  Pn  to  pre-Pn 
amplitude  ratio  (signal-to-noise)  cut-off  is  2:1  and  that  many  data  for  paths  crossing 
the  Iranian  Plateau  that  survive  this  criterion  for  lower  frequencies  (0.75-1.5  Hz)  are 
too  noisy  at  high  frequencies  (6.0-9.0  Hz).  This  suggests  that  attenuation  is  higher  in 
the  Iranian  lithosphere  relative  to  the  Kazahk  Platform.  For  the  lowest  frequency 
band,  0.75-1.5  Hz,  paths  crossing  the  Iranian  Plateau  generally  show  weak  Lg  (Pg/Lg 
>  1.0),  whereas  paths  crossing  the  Kazahk  Platform  (with  similar  distances  to  ABKT) 
show  strong  Lg  (Pg/Lg  <  0.5). 

The  complicated  regional  phase  behavior  of  this  region  makes  it  well  suited  for 
evaluating  regional  phase  measurement  techniques  because  a  wide  range  of  Pg/Lg 
ratios  is  observed.  However,  for  monitoring  purposes  azimuthal  variations  in  regional 
phase  amplitude  ratios  complicate  the  use  of  a  single  distance  correction.  Investigation 
of  the  azimuthal  dependence  of  regional  phase  propagation  at  this  station  and  distance 
corrections  will  be  presented  in  a  future  study. 

Time  Domain  Measurements  -  Methodology 


Time  domain  measurements  of  phase  amplitudes  were  made  by  narrow  band 
filtering  the  waveforms  using  a  4-pole,  2-pass  Butterworth  filter  (bands:  0.75-1.5;  1.5- 
3.0;  3.0-6.0;  6.0-9.0  Hz).  Pn,  Pg  and  Lg  phases  were  isolated  as  described  above.  For 
each  windowed,  de-meaned  phase,  s(t),  we  measured  the:  absolute  mean  (LI  norm): 

abs.  mean  =  / l(r)  ;  where  / l(r)  =  L (r)l  ;  (1) 


rms  (root-mean-square  or  L2  norm): 


rms  = 


(2) 
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ABKT  Pg/Lg 


Amplitude  Ratio 


Pg/Lg  <  0.5,  Lg  efficient 


Time  domain:  Rms  ♦ - 0.5  <=  Pg/Lg  <=  1.0,  Lg  intermediate 


3.0-6.0  Hz  6.0-9.0  Hz 


Figure  3.  Map  of  Lg  propagation  efficiency  for  station  ABKT.  Pg/Lg  ratios  plotted 
along  the  path  are  color  coded  (key  upper  right). 


envelope  mean: 


env.  mean  -  e  it  ) 


where  e(r)  = 


5(r)2  +  //[^(r)]2  ; 


(3) 


and  //[^(r)]  is  the  Hilbert  transform;  the  absolute  maximum,  abs.  max.  =  max[ll(t)]; 
and  envelope  maximum,  env.  max.  =  max[e(t)].  The  absolute  mean,  rms  and  envelope 
mean  are  measures  of  the  average  energy  in  the  window,  whereas  the  absolute  max¬ 
imum  and  envelope  maximum  are  measures  of  the  peak  energy.  Time  domain  meas¬ 
urements  for  each  phase  are  made  for  each  pass  band  separately. 


Time  Domain  Measurements  -  Results 


The  waveform  data  at  ABKT  were  recorded  as  digital  counts  proportional  to 
ground  velocity.  After  deconvolution  of  the  instrument  response  and  integration  the 
trace  represents  true  ground  displacement.  Comparisons  of  the  Pn/Lg  amplitude  ratios 
measured  in  the  time  domain  from  raw  velocity  and  instrument  deconvolved  displace¬ 
ment  seismograms  are  shown  in  Figure  4.  Rms  amplitudes  were  used  and  the  ratio  of 
Pn  to  the  pre-Pn  noise  was  required  to  be  2.0  or  greater.  The  comparison  is  shown  for 
the  same  four  frequency  bands  as  Figure  3.  The  individual  log^oCPn  /  Lg]  ratios  are 
plotted  as  circles  and  are  fit  to  a  straight  line.  Also  shown  in  each  panel  are  the 
number  of  points,  linear  correlation,  standard  error  about  the  best-fit  straight  line  and 
the  slope  of  the  fit.  The  raw  velocity  and  instrument  deconvolved  displacement  give 
very  similar  Pn/Lg  amplitude  ratios.  The  slopes  and  linear  correlations  are  all  close  to 
1.0.  As  noted  above,  the  number  of  data  surviving  the  signal-to-noise  criterion  rapidly 
decreases  as  the  frequency  pass  band  is  increased.  Also  note  that  the  lower  Pn/Lg 
ratios  (logio[Pn/Lg]<0.0)  are  preferentially  discarded  due  to  the  signal-to-noise  criteria. 

Comparison  of  the  various  time  domain  measurements  for  Pn/Lg  and  Pg/Lg  are 
complied  in  Tables  1  and  2,  respectively.  The  various  time  domain  measurements  give 
fairly  consistent  results.  The  absolute  and  envelope  measurements,  both  mean  and 
maximum,  are  particularly  strongly  related.  The  rms  (energy  measure)  and  the  abso¬ 
lute  maximum  (peak  measure)  are  not  correlated  as  well  as  the  other  pairs  shown.  The 
standard  error  generally  increases  and  the  linear  correlation  decreases  as  frequency 
increases. 

Frequency  Domain  Measurements  -  Methodology 
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Figure  4.  Comparison  of  Pn/Lg  amplitude  ratios  (rms  amplitude  in  the  time-domain) 
for  raw  velocity  and  instrument  deconvolved  displacement  for  the  four  pass  bands. 

The  pre-Pn  noise  to  Pn  amplitude  ratio  was  required  to  be  greater  than  or  equal  to  2.0. 

The  data  (circles)  were  fit  to  the  straight  (solid)  line.  1-a  errors  in  the  fit  are  shown  as 
dotted  lines.  The  pass  band,  number  of  data,  n,  linear  correlation,  r,  standard  error 
about  the  best-fitting  straight  line,  sd^and  the  slope  of  the  fit,  k,  are  given. 


Table  1. 


Pn/Lg  Amplitude  Ratios  -  Comparison 


Frequency  band 

6.0-9.0  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

109 

0.997 

0.019 

0.997 

abs.  mean  vs.  rms 

109 

0.977 

0.052 

1.032 

abs.  mean  vs.  env.  mean 

109 

1.000 

0.002 

1.001 

rms  vs.  env.  mean 

109 

0.977 

0.049 

0.926 

abs.  max.  vs.  env.  max. 

109 

0.998 

0.019 

0.996 

rms  vs.  abs.  max. 

109 

0.910 

0.127 

1.150 

rms  vs.  fireq.  log 

102 

0.717 

0.222 

0.823 

rms  vs.  freq.  lin. 

119 

0.789 

0.204 

0.990 

Frequency  band 

3.0-6.0  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

162 

0.978 

0.047 

0.973 

abs.  mean  vs.  rms 

161 

0.990 

0.032 

1.006 

abs.  mean  vs.  env.  mean 

161 

1.000 

0.002 

0.999 

rms  vs.  env.  mean 

161 

0.990 

0.031 

0.974 

abs.  max.  vs.  env,  max. 

161 

0.997 

0.021 

0.984 

rms  vs.  abs.  max. 

161 

0.931 

0.095 

1.060 

rms  vs.  freq.  log 

164 

0.889 

0.101 

0.864 

rms  vs.  freq.  lin. 

182 

0.955 

0.065 

0.927 

Frequency  band 

1.5-3.0  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

245 

0.991 

0.042 

0.937 

abs.  mean  vs.  rms 

245 

0.997 

0.026 

0.984 

abs.  mean  vs.  env.  mean 

245 

1.000 

0.001 

1.000 

rms  vs.  env.  mean 

245 

0.997 

0.026 

1.010 

abs.  max.  vs.  env.  max. 

245 

0.999 

0.013 

0.992 

rms  vs.  abs.  max. 

245 

0.967 

0.081 

0.955 

rms  vs.  freq.  log 

251 

0.953 

0.091 

0.890 

rms  vs.  freq.  lin. 

254 

0.979 

0.065 

0.984 

Frequency  band 

0.75-1.5  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

251 

0.997 

0.036 

0.974 

abs.  mean  vs.  rms 

251 

0.998 

0.028 

0.9% 

abs.  mean  vs.  env.  mean 

251 

1.000 

0.001 

1.000 

rms  vs.  env.  mean 

251 

0.998 

0.028 

1.001 

abs.  max.  vs.  env.  max. 

251 

1.000 

0.012 

1.000 

rms  vs.  abs.  max. 

251 

0.985 

0.079 

1.009 

rms  vs.  freq.  log 

256 

0.984 

0.080 

0.991 

rms  vs.  freq.  lin. 

254 

0.992 

0.056 

0.997 
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Table  2 


Pg/Lg  Amplitude  Ratios  -  Comparison 


Frequency  band  6.0-9.0  Hz 


Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

109 

0.997 

0.016 

0.988 

abs.  mean  vs.  rms 

109 

0.972 

0.046 

1.021 

abs,  mean  vs.  env.  mean 

109 

1.000 

0.002 

1.000 

rms  vs.  env.  mean 

109 

0.971 

0.044 

0.923 

abs.  max.  vs.  env.  max. 

109 

0.997 

0.020 

0.989 

rms  vs.  abs.  max. 

109 

0.875 

0.120 

1.127 

rms  vs,  freq.  log 

102 

0.625 

0.173 

0.689 

rms  vs.  freq.  lin. 

119 

0.714 

0.171 

0.888 

Frequency  band 

3.0-6.0  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

162 

0.971 

0.038 

0.975 

abs.  mean  vs.  rms 

161 

0.982 

0.030 

1.035 

abs,  mean  vs.  env.  mean 

161 

1.000 

0.002 

0.999 

rms  vs.  env.  mean 

161 

0.982 

0.029 

0.931 

abs.  max.  vs.  env.  max. 

161 

0.993 

0.024 

0.984 

rms  vs.  abs.  max. 

161 

0.909 

0.083 

1.134 

rms  vs.  freq.  log 

164 

0.832 

0.088 

0.810 

rms  vs.  freq.  lin. 

182 

0.892 

0.071 

0.871 

Frequency  band 

1. 5-3.0  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

245 

0.987 

0.035 

0.925 

abs.  mean  vs.  rms 

245 

0.995 

0.022 

1.014 

abs.  mean  vs.  env.  mean 

245 

1.000 

0.001 

0.999 

rms  vs.  env.  mean 

245 

0.995 

0.022 

0.977 

abs.  max.  vs.  env.  max. 

245 

0.999 

0.013 

0.988 

rms  vs,  abs.  max. 

245 

0.952 

0.075 

1.019 

rms  vs.  freq.  log 

251 

0.914 

0.084 

0.837 

rms  vs.  freq.  lin. 

254 

0.956 

0.065 

0.941 

Frequency  band 

0.75-1.5  Hz 

Measures 

n 

r 

sd 

k 

raw  vel.  vs.  deconed  dis. 

251 

0.994 

0.029 

0.977 

abs.  mean  vs.  rms 

251 

0.996 

0.026 

0.995 

abs.  mean  vs.  env.  mean 

251 

1.000 

0.001 

1.000 

rms  vs.  env.  mean 

251 

0.996 

0.026 

0.996 

abs.  max.  vs.  env.  max. 

251 

0.999 

0.012 

1.001 

rms  vs.  abs.  max. 

251 

0.961 

0.081 

1.009 

rms  vs.  freq.  log 

256 

0.968 

0.073 

1.007 

rms  vs.  freq.  lin. 

254 

0.980 

0.055 

0.996 
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Frequency  domain  measurements  were  made  by  windowing  the  noise  and  each 
phase  as  described  above,  de-meaning  and  tapering  the  trace,  then  computing  the  fre¬ 
quency  spectrum  for  a  broad  frequency  range  (0.03-9.99  Hz)  by  fast-Fourier  transform 
(fft).  The  frequency  spectrum  was  re-sampled  and  stored.  Then  the  logjo-mean  spec¬ 
tral  amplitude  was  computed  for  the  same  frequency  bands  considered  above  (0.75-1.5; 
1. 5-3.0;  3.0-6.0;  6.0-9.0  Hz).  Only  frequencies  within  the  pass  band  were  used  to 
compute  the  log^Q-niean  spectral  amplitude  of  each  phase.  An  immediate  advantage  of 
frequency  domain  techniques  is  that  the  broadband  spectra  of  each  phase  need  only  be 
calculated  once  and  the  operation  of  computing  the  log-mean  spectral  amplitude  is  fast. 

Comparing  Time  and  Frequency  Domain  Measurements 

Comparison  of  the  time  and  frequency  domain  Pg/Lg  phase  amplitude  ratios  is 
shown  in  Figure  5.  Statistics  for  comparison  of  time  and  frequency  domain  Pn/Lg  and 
Pg/Lg  amplitude  ratios  are  compiled  in  Tables  1  and  2,  respectively.  Time  domain 
measurements  were  made  on  instrument  deconvolved  displacement  data  using  the  rms 
measurement.  For  the  lowest  frequency  band  (0.75-1.5  Hz),  the  measurements  com¬ 
pare  quite  favorably,  although  they  are  more  scattered  than  most  of  the  time  domain 
comparisons.  However,  as  the  pass  band  increases  the  correlation  decreases,  the  slope 
deviates  from  1.0  and  the  data  become  more  scattered.  Scatter  is  partly  due  to  signal- 
to-noise.  We  found  that  increasing  the  signal-to-noise  criteria  reduced  the  standard 
errors,  but  the  linear  correlation  and  slope  were  largely  unaffected  and  the  number  of 
data  were  strongly  reduced.  Similar  decreases  in  correlation  and  slope  were  found  for 
Pn/Lg. 

The  trend  of  the  slopes  in  Figure  5  indicates  that  the  time  domain  measurements 
of  Pg/Lg  are  consistently  larger  than  the  frequency  domain  measurements  for  the 
higher  pass  bands.  Insight  into  possible  causes  of  this  disagreement  can  be  seen  by 
plotting  the  frequency  spectra  of  the  broadband  and  narrow  band  filtered  data,  shown 
in  Figure  6.  This  is  for  the  seismogram  shown  in  Figure  2.  The  noise  and  Pn,  Pg  and 
Lg  phases  were  windowed  as  described  above.  Then  the  broadband  (0.5-9.99  Hz)  and 
narrow  band  filtered  (0.75-1.5,  1. 5-3.0,  3.0-6.0  and  6.0-9.0  Hz,  2-,  4-  and  8-pole,  2- 
pass  Butterworth  filter)  spectra  were  computed  by  fft.  Within  the  0.75-1.5,  1. 5-3.0  and 
3. 0-6.0  Hz  pass  bands  the  spectra  of  the  filtered  data  compare  quite  well  with  the 
broadband  spectra.  However,  significant  energy  from  outside  the  pass  band,  particu¬ 
larly  at  frequencies  below  the  pass  band,  contributes  to  the  filtered  spectra.  The  time 
domain  filtering  is  equivalent  to  multiplication  in  the  frequency  domain  of  a  symmetric 
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Figure  5.  Comparison  of  Pg/Lg  amplitude  ratios  measured  in  the  time  (instrument 
deconvolved  displacement,  rms)  and  frequency  domains,  similar  to  Figure  4.  The  pre- 
Pn  noise  to  Pn  amplimde  ratio  was  required  to  be  greater  than  or  equal  to  2.0. 
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Figure  6.  Comparison  of  Pn,  Pg  and  Lg  spectra  for  broadband  (solid)  and  narrow  band 
filtered  data  (dashed)  for  the  waveform  shown  in  Figure  2.  The  noise  spectrum  is 
shown  as  the  dotted  black  line.  Spectra  are  plotted  as  logjQ  amplitude.  Vertical  lines 
indicate  the  filter  pass  band  frequencies. 
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Butterworth  filter  with  the  regional  phase  spectrum,  which  drops  off  rapidly  with 
increasing  frequency.  The  resulting  spectrum  is  skewed  to  lower  frequencies.  Increas¬ 
ing  the  number  of  poles  strongly  reduces  the  energy  from  outside  the  pass  band,  while 
maintaining  excellent  agreement  to  the  broadband  spectrum  within  the  pass  band.  We 
found  that  the  relationship  of  the  time  domain  measurements  to  the  frequency  domain 
measurements  was  generally  similar  to  that  shown  in  Figure  5,  regardless  of  how  many 
poles  were  used  for  the  time  domain  measurements.  This  suggests  that  energy  from 
outside  the  pass  band  is  not  the  primary  cause  of  the  disagreement  between  the  time 
and  frequency  domain  measurements. 

For  the  highest  frequency  band  (6.0-9.0  Hz),  note  that  the  broadband  spectra  in 
Figure  6  begin  to  increase  around  9.0  Hz.  This  results  from  deconvolution  of  the 
instrument  response.  The  anti-aliasing  filter  stans  rolling  off  at  around  6.0  Hz.  When 
the  observed  spectra  are  divided  by  this  response,  the  resulting  spectra  are  greatly 
amplified.  Clearly  there  are  significant  differences  between  the  broadband  and  filtered 
spectra  for  this  highest  pass  band.  The  noise  spectrum  is  also  shown  in  each  panel  of 
Figure  6  as  the  dotted  black  line.  For  the  Lg  phase,  noise  and  signal  spectra  merge  for 
the  low  end  of  this  highest  pass  band.  The  signal-to-noise  criteria  is  applied  to  the  Pn 
and  pre-Pn  time  windows.  For  the  case  shown,  the  Lg/pre-Pn  amplitude  ratio  would 
probably  fail  to  be  greater  than  2.0,  although  Pn  and  Pg  would  probably  pass.  Pn,  Pg, 
Sn  and  Lg  signal-to-noise  ratio  attributes  for  this  station  will  be  discussed  further  in  a 
future  paper. 

We  believe  that  the  main  cause  of  disagreement  between  the  time  and  frequency 
domain  measurements  in  Figure  5  is  due  to  log  averaging  of  the  spectra.  The  spectra 
for  all  phases  fall  off  very  rapidly  within  the  band  1.0-10  Hz.  Averaging  the  logjo  of 
the  spectrum  does  not  weight  the  amplitudes  associated  with  frequencies  near  the  low 
end  of  the  pass  band  as  heavily  as  does  linear  averaging  of  the  spectra.  Time  domain 
amplitude  measurements  weight  each  frequency  equally  and  are  thus  biased  toward  the 
relatively  larger  amplitude  energy  at  the  low  end  of  the  pass  band.  To  illustrate  this 
point,  we  computed  linear  averages  of  the  specn-al  amplitude  of  each  phase  and  formed 
the  amplitude  ratios  Pn/Lg  and  Pg/Lg.  These  measurements  were  well  correlated  with 
the  time  domain  (rms)  measurements  and  the  slopes  of  the  scatter  plots  were  closer  to 
1.0.  Results  are  compiled  in  Tables  1  and  2.  Linear  and  log  averaging  of  a  spectrum 
should  yield  converging  results  as  the  pass  band  is  narrowed.  Remaining  disagreement 
between  the  time  and  linearly  averaged  frequency  domain  measurements  for  the 
highest  pass  band  is  probably  due  to  poor  Lg  signal-to-noise  as  described  above. 


Discussion 


In  this  study  we  present  comparisons  of  regional  phase  amplitude  ratio  measure¬ 
ment  techniques.  Regional  phase  ratios  (Pn/Lg  and  Pg/Lg)  for  the  area  studied  show 
large  variation,  and  thus  results  for  this  region  should  be  applicable  to  other  regions. 
Various  time  domain  measurements,  particularly  energy  measures,  agree  very  well  with 
each  other  (linear  correlations  between  0.97-0.99).  Absolute  and  envelope  mean  and 
maxima  agree  nearly  perfectly.  Energy  and  peak  measures  do  not  compare  as  favor¬ 
ably  with  each  other,  but  they  are  still  rather  well  correlated  (linear  correlations 
between  0.88-0.99).  For  some  cases,  peak  measures  in  the  time  domain  may  discrim¬ 
inate  better  than  energy  measures.  This  could  be  due  to  the  impulsive  nature  of  P- 
waves  from  explosions.  The  population  studied  in  this  paper  are  all  earthquakes,  thus 
investigation  of  discriminant  performance  was  not  possible. 

Time  domain  (rms)  and  frequency  domain  measurements  compared  favorably  for 
the  lowest  frequency  band,  but  the  linear  correlation  decreased  as  the  pass  band 
increased.  More  importantly,  however,  is  that  the  time  domain  Pn/Lg  and  Pg/Lg 
amplitude  ratios  were  consistently  higher  than  the  log-averaged  frequency  domain 
ratios  and  that  the  disagreement  increased  as  the  pass  band  increased.  A  major  disad¬ 
vantage  of  time  domain  measurements  is  revealed  in  Figure  6.  Time  domain  filtering 
does  not  offer  much  control  over  the  frequencies  which  contribute  to  the  measurement. 
Clearly  significant  energy  from  frequencies  below  the  pass  band  contributes  to  the  time 
domain  measurements,  particularly  for  the  2-pole,  2-pass  filter.  For  the  highest  fre¬ 
quency  band  (6.0-9.0  Hz)  significant  differences  between  the  broadband  and  filtered 
spectra  (Figure  6)  results  from  deconvolution  of  the  instrument  response.  This  prob¬ 
ably  contributes  to  the  disagreement  seen  for  this  pass  band  in  Figure  5  and  Tables  1 
and  2. 

Linear  averaging  of  the  spectral  amplitude  brings  the  time  and  frequency  domain 
measurements  into  much  better  agreement.  Log-averaging  provides  a  more  robust  esti¬ 
mate  of  the  energy  within  a  given  pass  band  because  the  regional  phase  spectra  vary 
by  several  orders  of  magnitude  within  the  band  0.5-10.0  Hz.  Time  domain  measure¬ 
ments  weight  all  frequencies  within  the  band  equally  and  are  thus  biased  by  the  large 
spectral  amplitudes  associated  with  the  low  end  of  the  pass  band.  There  is  no  easy  way 
to  make  a  time  domain  measurement  that  is  not  affected  by  this  bias.  Furthermore, 
because  the  broadband  spectra  need  only  be  calculated  and  stored  once,  the  frequency 
domain  measurements  are  more  efficient.  Averages  of  the  spectral  amplitude  within 
any  band  can  be  easily  computed  from  the  broadband  spectra,  instead  of  having  to 
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read  the  data,  filter  and  then  compute  the  amplitude.  Thus,  we  believe  frequency 
domain  measurements  of  regional  phase  amplitudes  and  amplitude  ratios  are  superior 
to  time  domain  measurements. 
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Appendix 

WMQ  Discriminant  Measures  and  Path  Dependence 
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Figure  Ala.  Correlations  between  path  properties  for  the  WMQ  observations. 
Relationships  between  path  length  (distance)  and  surface  topography  parameters. 
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Figure  Alb.  Correlations  between  path  properties  for  the  WMQ  observations 
Relationships  between  mean  path  elevation  and  and  other  topographic  parameters. 
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Figure  Ale.  Correlations  between  path  properties  for  the  WMQ  observations. 
Relationships  between  mean  path  roughness  and  other  topographic  parameters. 
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Figure  Aid.  Correlations  between  path  properties  for  the  WMQ  observations. 
Relationships  between  path  length  (distance)  and  crustal  waveguide  parameters. 
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Figure  Ale.  Correlations  between  path  properties  for  the  WMQ  observations. 
Relationships  between  mean  Moho  thickness  and  surface  topography  parameters. 
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Figure  Alf.  Correlations  between  path  properties  for  the  WMQ  observations. 
Relationships  between  mean  sediment  thickness  and  surface  topography  parameters. 
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Figure  A2a.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  path  length.  CC  denotes  the  linear  conelation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A2b.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path.  CC  denotes  the  linear  correlation  coefficient  SD 
denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding 
to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2c.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path  after  correcting  for  the  distance  dependence  in 
Figme  A2a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A2d.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  stuface  roughness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 


75 


LOG[Pg/Lg]  LOG[Pg/Lg] 


LOG[Pg/Lg]  LOG[Pg/Lg] 

-2.5  -1.5  -0.5  0.5  1.5  2.5  -2.5  -1.5  -0.5  0.5 


0.3  0.8  1.3  1.8  0.3  0.8  1.3  1.8 


ROUGHNESS  (KM)  ROUGHNESS  (KM) 


0.3  0.8  1.3  1.8  0.3  0.8  1.3  1.8 


ROUGHNESS  (KM)  ROUGHNESS  (KM) 


Figure  A2e.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  Ala.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A2f.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  path  elevation.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2g.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  surface  roughness.  CC  denotes  the  linear 
ccttrelation  coefficient  SD  denotes  the  standard  deviation  of  Ae  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2h.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path.  CC  denotes  the  hn^  correlanon 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  imes 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2i.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A2a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A2j.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path.  CC  denotes  the  linear  correlation  coefficient. 
SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2k.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path  after  correcting  for  the  distance  dependence  in 
Figure  A2a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A21.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  il  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2m.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path  after  correcting  for  the  distance  dependence 
in  Figure  A2a.  CC  denotes  the  linear  correlation  coefficient  SD  denotes  the  standard 
deviation  of  the  linear  regression,  ivith  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A2n.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  minimum  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2o.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  maximum  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2p.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A2q.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path  after  correcting  for  tiie  distance 
dependence  in  Figure  A2a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  die 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 


88 


-2.5  -1.5  -0.5  0.5  1.5  2.5  -2.5  -1.5  -0.5  0.5 

LOG[Pg/Lg]  LOG[Pg/Lg] 


MAX.  SEDIMENTS  (KM)  MAX.  SEDIMENTS  (KM) 


MAX.  SEDIMENTS  (KM)  MAX.  SEDIMENTS  (KM) 


Figure  A2r.  Frequency  dependent  ratios  of  Pg/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  maximum  sediment  thickness  on  each  path.  CC  denotes  the  linear  correlahon 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  hnes 
corresponding  to  il  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3a.  Irequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  paA  length.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  +1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A3b.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path.  CC  denotes  the  linear  correlation  coefficient  SD 
denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding 
to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3c.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
nraction  of  mean  elevation  of  each  path  after  correcting  for  the  distance  dependence  in 
Figure  A3a.  CC  denotes  the  linear  correlation  coefficient,  SD  denotes  the  standard 
^viation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A3d.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
functicMi  of  mean  surface  roughness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3e,  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
f^cnon  of  me^  surface  roughness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A3a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A3f.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  path  elevation.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3g  Frequency  ^pendent  ratios  of  Pn/Lg  phases  recorded  at  WMO  plotted  as  a 
function  of  the  prc^uct  o^ath  length  times  mean  surface  roughness.  CC  denotes  the  linear 
coeffiaenL  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  hnes  corrcspondmg  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3h.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 


97 


LOG[Pn/Lg]  LOG[Pn/Lg] 


LOG[Pn/Lg]  LOG[Pn/Lg] 


WMQ  1.5-3  Hz 
Corrected 

CC=-0.2536 
SD=0.3697 
K=-1 3.2542 

— - - 1 

• 

• 

•  • 

• 

• 

-1.1  -0.6  -0.1  0.4  0.9  1.4  1.9  2.4 

SKEWNESS 


CM 


CM 


WMQ  3-6  Hz 

WMQ  6-9  Hz 

Corrected  k=-2.0068 

Corrected  k=-0.4660 

-  —• - —  ^  —  - - —  — • - •—  «-  __ 

•  *  •  *. - s^- 

\  ^  •  •  **!••  ••  ^ 

_ - - ^ -  ^ - - - 

•  1  • 

• 

•  •  •  • 

H - ! - 1 - 1 - \ - i - 1 - 

• 

- 1 - 1 - 1 - 1 - 1 - 1 - r 

-1.1  -0.6  -0.1  0.4  0.9  1.4  1.9  2.4 

SKEWNESS 


-1.1  -0.6  -0.1  0.4  0.9  1.4  1.9  2.4 

SKEWNESS 


Figure  A3i.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A3a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A3j.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path.  CC  denotes  the  linear  correlation  coefficient. 
SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  il  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3k.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path  after  correcting  for  the  distance  dependence  in 
Fig^  A3a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  +1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A31.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3m.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
fiinction  of  mean  crustal  thickness  of  each  path  after  correcting  for  the  distance  dependence 
m  Figure  A3a.  CC  denotes  the  linear  correlation  coefficient  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression.  ’ 
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Figure  A3n.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  ^  a 
function  of  minimum  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3o.  Fr^uency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMO  plotted  as  a 
to^on  of  m^um  crustal  thiclmess  on  each  path.  CC  denotes  the  lin^^  collation 
oeffiaenL  SD  denot^  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3p.  Requency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A3q.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A3a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A3r.  Frequency  dependent  ratios  of  Pn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  maximum  sediment  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4a.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  path  length.  CC  denotes  die  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A4b.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path.  CC  denotes  the  linear  correlation  coefficient  SD 
denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding 
to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4c.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plo^  as  a 
function  of  mean  elevation  of  each  path  after  correcting  for  the  jis^ce  de^ndence  m 
Figure  A4a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  stand^ 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A4d.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4e.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A4a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A4f.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  path  elevation.  CC  denotes  the  linear 
correlation  coefficienL  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4g.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  surface  roughness.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4h.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4i  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
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and  K  denotes  the  slope  of  the  regression.  ^  ^ 
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Figure  A4j.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  ^  a 
function  of  mean  surface  slope  of  each  path.  CC  denotes  the  linear  correlation  coefficient. 
SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  +1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4k.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path  after  correcting  for  the  distance  dependence  in 
Figure  A4a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A41.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4m.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
fonction  of  mean  crustal  thickness  of  each  path  after  correcting  for  the  distance  dependence 
in  Figure  A4a.  CC  denotes  the  linear  correlation  coefficient  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A4n.  I^uency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  minimum  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4o.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  inaxiniuni  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4p.  Rtequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A4q.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
^notion  of  mean  segment  thickness  of  each  path  after  correcting  for  the  distance 
dependence  m  Figure  A4a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
stantod  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  +1  SD 
and  K  denotes  the  slope  of  the  regression.  b  -  ^  , 
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Figure  A4r.  Frequency  dependent  ratios  of  Pg/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  maximum  sediment  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5a.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  path  length.  CC  denotes  the  linear  coirelation  coefficient.  SD  denotes  die 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A5b.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path.  CC  denotes  the  linear  correlation  coefficient  SD 
denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding 
to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5c.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 

after  correcting  for  the  distance  dependence  in 
^  f  S-  correlation  coefficient.  SD  denotes  the  standard 

^viation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  +1  SD  and  K 
denotes  the  slope  of  the  regression.  &  -  d. 
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Figure  A5d.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5e.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A5a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  tfie 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A5f.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  path  elevation.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5g.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  pnxiuct  of  path  length  times  mean  surface  roughness.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5h.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5i.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path  after  correcting  for  the  distance 
dependence  in  Figure  A5a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A5j.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  ^  a 
function  of  mean  surface  slope  of  each  path.  CC  denotes  the  linear  correlation  coefficient. 
SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5k.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
firaction  of  mean  surface  slope  of  each  path  after  correcting  for  the  distance  dependence  in 
Fig^  A5a.  CC  denotes  the  linear  conelation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A51.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  +1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5m.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  crustal  thickness  of  each  path  after  correcting  for  the  distance  dependence 
in  Figure  A5a.  CC  denotes  the  linear  correlation  coefficient  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A5n.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  minimuni  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5o.  Frequency  dependent  ratios  of  Pn/Sn  phases  reconied  at  WMQ  plotted  as  a 
function  of  maximum  crustal  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient.  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5p.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  il  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A5q.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  sediment  thickness  of  each  path  after  correcting  for  the  distance 
dependence  m  Figure  A5a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A5r.  Frequency  dependent  ratios  of  Pn/Sn  phases  recorded  at  WMQ  plotted  as  a 
function  of  maximum  sediment  thickness  on  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A6a.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
fimetion  of  path  length.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD, 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A6b.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path.  CC  denotes  the  linear  correlation  coefficient  SD 
denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding 
to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 


145 


LOG[Sn/Lg]  LOG[Sn/Lg] 


LOG[Sn/Lg]  LOG[Sn/Lg] 


0.5  1.0  1.5  2.0  2.5  3.0  3.5  0.5  1.0  1.5  2.0  2.5  3.0  3.5 


MEAN  ELEVATION  (KM)  MEAN  ELEVATION  (KM) 


0.5  1.0  1.5  2.0  2.5  3.0  3.5  0.5  1.0  1.5  2.0  2.5  3.0  3.5 


MEAN  ELEVATION  (KM)  MEAN  ELEVATION  (KM) 


Figure  A6c.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  elevation  of  each  path  after  correcting  for  the  distance  dependence  in 
Fig^  A6a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the  standard 
deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD,  and  K 
denotes  the  slope  of  the  regression. 
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Figure  A6d.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  roughness  of  each  path.  CC  denotes  the  ^ear  correlanon 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  Imes 
corresponding  to  il  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A6e.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
Action  of  mean  surface  roughness  of  each  path  after  correcting  for  the  distance 
dependence  m  Figure  A6a.  CC  denotes  the  linear  conelation  coefficient  SD  denotefke 

deviation  of  the  hnear  regression,  with  the  dotted  lines  coiresponding  to  +1  SD 
and  K  denotes  the  slope  of  the  regression.  ^  'S  -a  ciu. 
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Figure  A6f.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  path  elevation.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A6g.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  the  product  of  path  length  times  mean  surface  roughness.  CC  denotes  the  linear 
correlation  coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the 
dotted  lines  corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figiue  A6h.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  skewness  of  each  path.  CC  denotes  the  linear  correlation 
coefficient  SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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Figure  A6i  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  me^  surface  skewness  of  each  path  after  correcting  for  the  distance 
dependence  m  Figure  A6a.  CC  denotes  the  linear  correlation  coefficient.  SD  denotes  the 
stantod  deviation  of  the  linear  regression,  with  the  dotted  lines  corresponding  to  ±1  SD 
and  K  denotes  the  slope  of  the  regression. 
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Figure  A6j.  Frequency  dependent  ratios  of  Sn/Lg  phases  recorded  at  WMQ  plotted  as  a 
function  of  mean  surface  slope  of  each  path.  CC  denotes  the  linear  correlation  coefficient. 
SD  denotes  the  standard  deviation  of  the  linear  regression,  with  the  dotted  lines 
corresponding  to  ±1  SD,  and  K  denotes  the  slope  of  the  regression. 
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